Kermions: quantization of fermions on Kerr space-time 
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We study a quantum fermion field on a background non-extremal Kerr black hole. We discuss 
the definition of the standard black hole quantum states (Boulware, Unruh and Hartle-Hawking), 
focussing particularly on the differences between fermionic and bosonic quantum field theory. Since 
all fermion modes (both particle and anti-particle) have positive norm, there is much greater flex- 
ibility in how quantum states are defined compared with the bosonic case. In particular, we are 
able to define a candidate 'Boulware'-like state, empty at both past and future null infinity; and 
a candidate 'Hartle-Hawking'-like equilibrium state, representing a thermal bath of fermions sur- 
rounding the black hole. Neither of these states have analogues for bosons on a non-extremal Kerr 
black hole and both have physically attractive regularity properties. We also define a number of 
other quantum states, numerically compute differences in expectation values of the fermion current 
and stress-energy tensor between two states, and discuss their physical properties. 
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I. INTRODUCTION 

In the absence of a definitive theory of quantum grav- 
ity, it is appropriate to attack the problem from a variety 
of directions. Quantum field theory in curved space-time 
treats the space-time geometry as a fixed, classical back- 
ground described by Einstein's field equations of general 
relativity. The behaviour of quantum matter fields on 
this background is then studied. This may be regarded 
as a first approximation to a full theory of quantum grav- 
ity (in which both the geometry and matter fields would 
be quantized). 

Central to the study of quantum fields on any partic- 
ular space-time background is the concept of a vacuum. 
For a free quantum field, the field is typically decom- 
posed into an orthonormal basis of positive and negative 
frequency field modes. In general the split into posi- 
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tive and negative frequency modes is not unique, par- 
ticularly when the background space-time does not pos- 
sess a globally-timc-like Killing vector. For a fixed split- 
ting of the quantum field into positive and negative fre- 
quency modes, the coefficients of the positive and neg- 
ative frequency modes are promoted to operators. The 
coefficients of the positive frequency modes become par- 
ticle annihilation operators and those of the negative 
frequency modes become particle creation operators. A 
'vacuum' state is defined as that state annihilated by the 
particle annihilation operators. The non-uniqueness of 
the splitting into positive and negative frequency modes 
therefore leads to a non-uniqueness of the definition of 
'vacuum'. For a general space-time, and for black hole 
space-times in particular, there may be several quan- 
tum states of physical interest which arise as 'vacuum' 
states from different ways of splitting the quantum field 
into positive and negative frequency modes. Even in 
Minkowski space, the concept of a 'vacuum' is observer- 
dependent, as demonstrated by the Unruh effect 

We now describe the main quantum states specifically 
on a Schwarzschild black hole background, since it is on 
this background where the states were originally defined 
and where their properties are better established [|| . 
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• The Unruh state Q models a spherically- 
symmetric, evaporating black hole formed by gravi- 
tational collapse. The Unruh state is empty at past 
null infinity, containing a quantum flux of thermal 
Hawking radiation emitted away to future null in- 
finity. While the Unruh state is irregular at the 'un- 
physical' past horizon, it is regular at the 'physical' 
future horizon. This state is clearly not invariant 
under the Schwarzschild symmetry of time-reversal, 
as the process of gravitational collapse itself is not 
time-reversal invariant. 

• The Hartle-Hawking state [f| represents a black 
hole in unstable thermal equilibrium with a bath 
of quantum radiation at the Hawking temperature. 
The Hartlc-Hawking state is particularly important 
in that it respects the symmetries of the underlying 
Schwarzschild space-time and is regular everywhere 
on and outside the event horizon. It is therefore the 
relevant state for black hole thermodynamics (see, 
for example, Q). Furthermore, physically, it is the 
state which is seen as empty by a freely-falling ob- 
server near the event horizon Q and, practically, 
this state is the easiest one to renormalize (see, 
for example, [H H, Q). We note that the equiv- 
alent of this state in Schwarzschild- AdS (anti-de 
Sitter) space-time is the one which is of relevance 
for black hole thermodynamics [l(| in that case and 
so for considering black holes in the context of the 
AdS/CFT (conformal field theory) correspondence 

SEMI- 

• The Boulware state [l4| models not a black hole 
but a (static and spherically-symmetric) cold star: 
it is divergent on the horizon (both future and 
past) and it is empty at radial infinity (both future 
and past). This state respects the symmetries of 
the Schwarzschild space-time, in particular, time- 
reversal symmetry. 

We note that, in Schwarzschild, the properties of the 
above states are the same independently of whether the 
quantized field is bosonic or fermionic [H, 0, [l4| . 

Our focus in this paper is the quantization of fermion 
fields on a non-extremal Kerr black hole background. 
The study of quantum fields propagating on a Kerr black 
hole has a long history, the discovery of 'quantum super- 
radiance' (the 'Unruh-Starobinskii' effect [HI, [l6|) pre- 
dating the famous Hawking radiation. However, apart 
from computations of the fermion Hawking flux from 
a Kerr black hole [l7l - [20j or on-the-brane emission of 
fermions from a higher-dimensional rotating black hole 
pH l22j , most of the work in the literature has focussed on 
bosonic quantum fields on Kerr. A key feature of classical 
bosonic fields on Kerr is super-radiance [23| , whereby an 
incoming wave can be reflected back to infinity with an 
amplitude greater than initially. In contrast, fermionic 
fields do not exhibit classical super-radiance [23| (we 
note, however, that a classical fermion field might not 



have a clearly well-defined physical meaning [24| , and use 
the term 'classical' to denote a field which is not quan- 
tized and satisfies a wave equation). Quantum super- 
radiance (the 'Unruh-Starobinskii' radiation) is nonethe- 
less present for fermions as well as bosons [HI, [l6| . This 
lack of classical super-radiance for fermion fields is one 
motivation for our investigation of the properties of quan- 
tum fermion fields on a Kerr black hole. 

Quantum scalar fields have received particular atten- 
tion. Notable is the theorem of Kay and Wald plj . 
proved for scalar fields, that there docs not exist a 
Hadamard state (that is, a state whose short-distance 
singularity structure is of the Hadamard form - see, for 
example, [111) on Kerr which is regular everywhere and 
preserves the symmetries of the space-time. This means, 
in particular, that there is no analogue of the 'Hartle- 
Hawking' state in the Schwarzschild space-time Q for 
scalar fields on Kerr. While there have been attempts in 
the literature to define a state for bosons which mimics at 
least some of the properties of the Hartle-Hawking state 
0, H3] , these states either do not represent an equilibrium 
state or fail to be regular almost everywhere [28|, [29| . In 
particular, the Frolov-Thorne state Q, constructed using 
the 7?-formalism, is regular only on the axis of rotation of 
the black hole |28| , and is ill-defined everywhere else even 
inside the speed-of-light surface (defined in Sec. Ill Al) . 
A solution is to place a mirror inside the speed-of-light 
surface, and then a regular equilibrium thermal state re- 
specting the symmetries of the space-time geometry in- 
side the mirror can be constructed [3(| ■ 

For both scalar [28[ and electromagnetic fields [29[ in 
Kerr a 'past-Boulware' state can be constructed, which 
is empty at past null infinity Z~ but not at future null 
infinity T + (see Fig. [T]), where it contains the 'quantum 
super-radiance'. Numerical computations of differences 
of expectation values in this state and the 'past-Unruh' 
state [28[ (which is empty at I~ , contains the Hawking 
radiation at I + and is the analogue for Kerr black holes 
of the Unruh state [|[ for Schwarzschild black holes) for 
electromagnetic fields can be found in (29|. The lack 
of an analogue in Kerr of the 'Hartle-Hawking' state in 
Schwarzschild for bosonic fields is linked to a similar lack 
of a true 'Boulware' state which is empty at both T~ and 

With such a consistent picture developed for both 
scalars and electromagnetic radiation, there may seem 
to be little merit in a detailed study of the quantum 
field theory of fermions on Kerr, which is perhaps why 
none has been attempted to date. However, we will 
show that quantum fermion fields are rather different to 
quantum bosonic fields on Kerr black holes. In partic- 
ular, the lack of classical super-radiance makes the de- 
velopment of canonical quantization rather simpler for 
fermions than for bosons. However, the differences are 
not simply technical, but deeper as well. We are able to 
define analogues of the 'Hartle-Hawking' [j| and 'Boul- 
ware' [l4| vacua which are closer approximations to the 
corresponding states on Schwarzschild space-time than 
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is possible for bosonic fields on Kerr. The new fermionic 
states that we define have divergences which can never- 
theless be understood physically: the 'Hartle-Hawking' 
state diverges on and outside the speed-of-light surface 
(in the region where an observer co-rotating with the 
event horizon must have a velocity greater than or equal 
to the speed of light) and the 'Boulware' state diverges 
in the ergosphere (the region where an observer cannot 
remain at rest with respect to infinity - see Sec. Ill A[) . 



The outline of this paper is as follows. In Sec. |TT] we 
review the salient features of the Kerr space-time and 
the classical mode solutions of the Dirac equation on this 
background. The canonical quantum theory of fcrmions 
on Kerr is developed in Sec. IIII1 where we focus in par- 
ticular on defining quantum states, firstly the uncon- 
troversial 'past-Boulware' and 'past-Unruh' states, and 
secondly we present candidate 'Boulware' and 'Hartle- 
Hawking' states. The properties of these states are in- 
vestigated in Sec. [TV] where we compute the differences 
in expectation values of the fermion number current and 
stress-energy tensor in two different states. The lack 
of a suitable renormalization procedure for fermions on 
Kerr (unlike that for Schwarzschild 0, [3l[ ) means that 
differences in expectation values between two states are 
all that are currently tractable. Our conclusions on the 
physical properties of the states we have constructed are 
summarized in Sec. [V] The implications of our results 
are discussed in Sec. I VII including their relevance to the 
Kerr-CFT correspondence [HJ (see also [33|, [34| for re- 
views) . 



II. SPIN-1/2 PARTICLES ON KERR 

A. Kerr geometry 

The Kerr metric in the usual Boyer-Lindquist co- 
ordinates (t, r, 9, ip) has the form 

ds 2 = -— [dt - a sin 2 6 dip] 2 + —dr 2 + £ d6 2 



sin r , 2 2 \ , ,,-,2 
H — — [{r 2 + a 2 ) dtp - adt\ 



FIG. 1. Part of the Carter-Penrose diagram for the complete where 
Kerr geometry, showing the future event horizon H + , past 
event horizon H~ , future null infinity X + and past null infin- A 
ity Region I corresponds to the space-time exterior to 
the event horizon and is the region on which we study the 
quantum fermion field. Region IV will be required in Sec. lIIII 
for defining some of our quantum states. A more complete 
Carter-Penrose diagram for the Kerr geometry can be found 
in [H. 
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S = r 2 + a 2 cos 2 i 



(2.1) 



(2.2) 



with M the mass of the black hole and J = aM its 
angular momentum. Here, and throughout this paper, we 
use units in which c = G = H = ks = 1- We employ the 
space-time signature ( — h H — h), which means that care 
has to be taken, particularly with the Dirac matrices (|B1[) 
and spin connection matrices (|B8|) . because many papers 
in the quantum field theory literature use the alternative 

signature (H ). 

The outer event horizon of the Kerr black hole is at 



M + y/M 2 - a 2 



and has Hawking temperature 



T H = 



H 



(2.3) 



(2.4) 



u 



In this paper we consider only non-extremal Kerr black 
holes, for which the outer event horizon has non-zero 
Hawking temperature and < a < M. Part of the 
Carter-Penrose diagram of the full non-extremal Kerr 
space-time is shown in Fig. [TJ 

The Kerr metric (|2.ip is stationary and axisymmetric, 
possessing two Killing vectors: 



d_ 



x 



d_ 

dtp' 



(2.5) 



The Killing vector £ is time-like near infinity, but be- 
comes null on the surface given by 



rs 



AI+ \JM 2 - a 2 i 



(2.6) 



namely the stationary limit surface. Inside the station- 
ary limit surface (the region between the stationary limit 
surface and the event horizon being the ergosphere) , the 
vector £ is space-like, indicating that, inside the ergo- 
sphere, observers cannot remain at rest relative to infin- 
ity. For a non-extremal black hole, the alternative Killing 
vector 



where 



H 



(2.7) 
(2.8) 
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FIG. 2. The cross-section of the stationary limit (red) and speed-of-light (blue) surfaces, for a < ao = Afy2 [>/2 — l] (left), 
a = ao (centre) and a > ao (right). In each case we have plotted cross-sections on a plane of fixed azimuthal angle cp. The axis 
of rotation of the black hole is a vertical line through the centre of each diagram, and the equatorial plane a horizontal line 
through the centre of each diagram. The black circle denotes the region inside the event horizon. 



is the angular velocity of the event horizon, is time-like 
sufficiently close to the horizon, becoming null on the 
event horizon (of which it is the generator). The Killing 
vector £ remains time-like outside the event horizon up to 
the speed-of-light surface (which we denote Sl), on which 
it becomes null. Physically, Sl is the surface outside 
which an observer can no longer have the same angular 
velocity as the event horizon. 

The surface Sl is distinct from the stationary limit 
surface and its location is given by the solution of a cu- 
bic equation for r in terms of 9, which can be found in 
the Appendix of [3(|. The smallest value of r on Sl 
arises in the equatorial plane 6 = ^ , while r — > oo on Sl 
as 9 — > 0, 7r and the axis of rotation is approached. In 

[30j . it is shown that for a < MJ2 \y/2 — l] the speed- 
of-light surface lies entirely outside the ergosphere; for 
M^2 [\/2 — l] < a < M part of Sl near the equato- 
rial plane lies inside the ergosphere. When a — ao = 
M J 2 [V2 — l] , the stationary limit surface touches the 
speed-of-light surface on the circle at r = 2M, 9 = ^ . For 
an extremal black hole a = M, the speed-of-light surface 
touches the event horizon in the equatorial plane. The 
location of the stationary limit surface and speed-of-light 
surface is shown in Fi g. [2] for the cases a < ao , a = ao 
and a > ao (see also |36[ for a recent discussion of the 
speed-of-light surface for Kerr). 



B. Formalism for fermions in curved space 

We consider massless fermions of spin- 1/2 propagating 
on the fixed Kerr geometry (12.11) . We use Dirac 4-spinors 
and our formalism follows [l5|, modulo some changes of 
sign due to our different convention for the space-time 



signature. We restrict our attention to massless fermions 
for simplicity. While the formalism developed in this 
and the following subsection is standard [l5|, H^-El , we 
explicitly give all our definitions so that our presentation 
is consistent and to avoid confusion. 

We begin with the Dirac equation for massless fermions 
on the Kerr space-time: 

7 ^V M * = 0, (2.9) 

where ^ is a Dirac 4-spinor. The Dirac matrices 7 M sat- 
isfy the anti-commutation relations 

7 Y+7Y = 2r, (2-10) 

where is the inverse metric. A suitable basis of 7^ 
matrices for the Kerr metric (|2.ip can be found in [l^, 137} 
and is reproduced for convenience in App. [B] Except in 
App. [5] throughout this paper the operators V M are the 
spinor covariant derivatives defined in terms of the spinor 
connection matrices r M as follows [TEl ]: 

V ^ = ^- T ^- (2 - U) 

The spinor connection matrices V ^ are defined in terms 
of covariant derivatives of the Dirac matrices 7^: 

+ r^ 7 K - r,7 M + 7T„ = 0, (2.12) 

where T^ K are the usual Christoffel symbols. A suitable 
choice of the spinor connection matrices for the Kerr 
metric can be found in App. [B] 

Massless fermion solutions to the Dirac equation (|2.9[) 
can be classified as "left-handed" or "right-handed" as 
follows. We first define a chirality matrix 7 s by 

7 5 = k^/l^'l 1 , (2-13) 
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where e^ u \ a is the Levi-Civita anti-symmetric symbol 
and i = \J— 1. The form of j 5 can be found in App. [5] 
Spinors are "left-handed" if they satisfy the equation 

MM 



(l -7°) f = 
and "right-handed" if they satisfy 
(1 +7 5 ) f = 0. 



(2.14) 



(2.15) 



If \I/ is a solution of the Dirac equation (12.91) . then 7 2 \£'* 
is also a solution of the Dirac equation |38[, where j 2 is 
a flat-space Dirac matrix given in App. [B] and the as- 
terix denotes complex conjugation. Furthermore, if 
is a left-handed spinor, then j 2 ^/* is right-handed, and 
vice- versa. 

The action giving rise to the field equation (|2.9|) is 



5 



d A x, 



-9 [*VV M * - (V„tt) 7 M *] (2-16) 



where the conjugate spinor 'J is given by 'J = ^a, with 
"ft the usual hermitian conjugate of "J considered as a 
matrix. The matrix a satisfies the conditions 



= ajV + 7 Mt a, 







(2.17) 



and a suitable choice of a is simply a = — 7 where 7 is a 
flat-space Dirac matrix defined in App. |Bl Note that this 
definition of the matrix a involves a minus sign relative 
to much of the literature, due to our metric conventions. 
The covariant derivative of the conjugate spinor W is 



(2.18) 



From the action (12.161) the classical stress-energy ten- 
sor is obtained 140, \4M : 



T ^ = 2 [^70, V„)* - (V (M W) 7,)*] , (2.19) 

where parentheses are used to denote symmetrization of 
indices. 

For any two spinor solutions of the Dirac equation 
(|2.9p . VPi and ^2, we define a conserved current J M fijl: 



J" = ^17^^ 



(2.20) 



An inner product between two solutions may be defined 
with respect to a constant t hypersurface St using the 
current component J*, as follows: 



(*i,^ 2 ) = / ^i^n^ 2 dS, 



S, 



(2.21) 



where n M is the unit outwards-pointing normal to St 



C. Solutions of the Dirac equation on Kerr 

The Dirac equation (12.91) is known to be separable on 
the Kerr geometry [38l l42j . Mode solutions take the form 
0,11113: 



'/A 

Lr] A 



(2.22) 



Spinors with L = +1 arc "left-handed" while those with 
L = — 1 are "right-handed" . The function T in (|2 .22[) is 
given by (37| 



T = 



A (r — iaL cos 9) sin 2 9 



(2.23) 



where we have corrected a sign error which appears in 
many places in the literature. The two-spinor 7/a is 



1 R A (r) 1 S A (9) 
2RA(r) 2 S A {9) 



(2.24) 



where A = {uj,£, m} is the set of quantum numbers for 
each spinor mode. Throughout this paper, the quanti- 
ties co, £, m and therefore ui = u> — miln are real; the 
quantities £ and m are half-integers. 

The radial and ang ular functions satisfy, respectively, 
the equations [l5l.ll9l [38j: 



d iKL 
dr A~~ 

d iKL 



li?A = A 2-Ra, 
2 i?A = A ii?A, 



where A is a separation constant (with A 
I = |, §, . . . when a = 0), 



(2.25) 



I for 



A" 



a 2 ) oj 



and 



d ( . . m 

— + ousm«- -— 
d9 V sm 

d / m 

— - awsin«- -— 
rifl V sm 9 



n, (2.26) 

i5a = A 2<Sa, 

2 5 A = -A 1 5 A . (2.27) 



It should be noted that the angular functions i/ 2 S\ are 
real but the radial functions \/ 2 Ra ar e complex. The ra- 
dial equations (|2.25[) depend explicitly on L. From (|2.25[) . 
under the mapping L — > — L the ordinary differential 
equations satisfied by the radial functions iRa and 2 Ra 
are interchanged. In our discussion below of particular 
mode solutions of the radial equations, we will be impos- 
ing boundary conditions on the radial functions which 
are valid for L = +1 only. The corresponding boundary 
conditions for L = — 1 can be found by swapping \Ra 
and 2 Ra- This should be borne in mind in later sections 
where physical quantities will depend on \Ra and 2-Ra- 
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Using the notation —A = {— uj,£,— m}, the following 
symmetries of the radial and angular functions will be 
useful for later calculations: 



iR- 



lR*A, 



iR- 



p* 



and 



lS-A — i2<?A) 



(2.28) 



(2.29) 



In (|2.29[) there is an ambiguity in an overall sign, which is 
irrelevant for the computation of physical quantities and 
can be chosen arbitrarily. The angular functions have an 
additional symmetry under 9 — > tt — 8: 

iSa (tt - 0) = ± 2 S A (0) 1 2 s a (k-o) = ±iS A (e). 

(2.30) 

We normalize the angular functions so that 



iS A (9) 2 de 



2 S A {9) 2 d6 = 1. 



(2.31) 



If i/ja is a solution of the Dirac equation (|2.9p . then so 
too is ip-A- However, we note that, despite the relations 
(|2.28l I2.29j) . ip-A is not equal to i/j a because of the com- 
plex function T (|2.23|) . If tpA is a solution of the Dirac 
equation with L = +1, then we can construct a corre- 
sponding solution with L = — 1 by changing L in (|2.22l 
I2.23P and in the radial equations (|2.25|) . 

The norm of the modes ()2.22j) is given by HV'aII = 
(V'AjV'a), where the inner product is given by (|2.21[) . It 
can be shown that all modes (|2.22|) have positive norm, 
regardless of the values of any of the quantum numbers. 
This is in contrast to the scalar case, where the sign of 
the norm of the modes depends on the frequency u) and 
azimuthal quantum number m (see App. [A"Tj) . 

One basis of mode solutions to the radial equations 
p.25p can be formed from the usual "in" and "up" radial 
functions (the expressions below are for the L = +1 case, 
the expressions in the L = — 1 case are found by making 
the transformation \Ra <-> 2 Ra)'- 



(iRa, 2 Ra) - 



(lR7,2R7) 



00 



(2.32) 



r* — > —00 
r\ —> 00 



(2.33) 



where uj = lu — rnQn and we have introduced the usual 
'tortoise' co-ordinate r*, defined by 



dr* r 2 + a 2 
~aV = A : 



(2.34) 



so that — > —00 at the event horizon and r* — > 00 as 
r — > 00. 

We also introduce an alternative basis, namely the 
"out" and "down" radial functions (as above, these ex- 
pressions are for the L = +1 case, swapping \Ra and 



2-Ra gives the expressions for the L = — 1 case): 



I pout po 



( pdown pdown\ 
U A > 2-Ka J 



'A 

Sur* /lout, 

, J± A t 



— > — 00 
— > 00 



(2.35) 



4< 



own iwr* 



— > —OO 
— > OO. 



(2.36) 



Unlike the scalar case (see (|A1|) in App.|X|, for fermions 
there are no particular subtleties in defining the "up" or 
"down" modes. This is because all the "up" and "down" 
modes have positive norm, independent of the sign of Q. 
This is our first indication that quantum field theory of 
fermions on Kerr may be more straightforward than that 
for bosonic fields. 

For any two solutions (iRa, 2 Ra) and ^i-Ra,2-Ra^ of 
the radial equations (|2.25[) . the quantities 

Wi = iR A 2 Ra - 2 Ra iRa, 

VK 2 = 1J R A1 i?A- 2 i?X 2 i? A , (2.37) 

can be shown to be independent of r. These two quan- 
tities can be used to derive a number of relationships 
between the constants in the functions (|2.32[ 12.331 12.351 
I2.36p . for example: 



1 



A1 



.in/up 
^A 



\A7\ 



l/up 



(2.38) 



with similar relations holding for "out / down" . From 
(|2.38p . we see that \A A \ < 1 for all modes, so that 
there is no classical super-radiance for fermions [23j (com- 
pare (|A5[) for the scalar case). To understand this lack 
of classical super-radiance for fermions, it is important 
to note that classical fermion fields do not satisfy the 
weak energy condition (the weak energy condition being 
that T^uU^y? > where it M is the four-velocity of any 
physical observer) [2^, |43| . Super-radiance for classical 
bosonic fields can be deduced from the area theorem be- 
cause these fields do satisfy the weak energy condition 
[23l . The fact that fermionic fields do not satisfy the 
weak energy condition means that they do not necessarily 
have to exhibit super-radiance, because the area theorem 
no longer holds. For the quantum field theory of bosonic 
fields on Kerr black holes, the existence of super-radiant 
modes causes many technical and conceptual difficulties 
0, l28l - [3"o1 | . While there is no classical super-radiance for 
fermions, it is still the case that the frequency of the 
modes as seen by an observer near infinity is u, while for 
an observer near the event horizon it is a), so subtleties 
remain. Despite the lack of classical super-radiance for 
fermions, we shall still use the terminology 'super-radiant 
modes' for those fermion modes for which ww < (which 
is the condition for super-radiance for scalar field modes, 
see Sec. lAip . 
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The "out" ([235)1 and "down" ([236]) radial functions 
can be compactly written in terms of the "in" (|2.32[) and 
"up" (|2.33[) radial functions as follows: 



nout 

1,2-n-A 

ndown a down r> 

1,2-Ka ^ A A 1,2-Ka 



/lout r>in i nout d u P 
A\ 1,2-rtA + -°A 1,2-Ka ' 



up 



B 



down 



pin 
1,2-flA 



(2.39) 



and the "in" and "up" radial functions can similarly be 
written in terms of the "out" and "down" radial func- 
tions. One important point for our later work is that the 
relations f|2 . 39[) only involve Ra, and not R A . This is in 
contrast to the situation for scalar fields, see App. [A] 

By inserting the appropriate radial functions into the 
two-spinor tja (|2.24[) we can construct basis spinor modes 

V>XS -0a\ ^A ut and V4° wn ( see SecEU. The "in" modes 
tp A correspond to unit flux incoming from past null in- 
finity T~ , part of which is scattered back to future null 
infinity T + and part passes down the future event hori- 
zon H + (see Fig. [1]). The "up" modes ip A p correspond to 
unit flux outgoing from the past event horizon TAT , part 
of which is scattered down the future event horizon % + 
and the rest travels out to I + . The "out" and "down" 
modes are the time reverse of the "in" and "up" modes: 
the "out" modes ip A ut correspond to unit flux outgoing 
at future null infinity X + , part of which has come from 
past null infinity T~ and part from the past event horizon 
7i~ . Similarly, the "down" modes ^ own correspond to 
unit flux going down the future event horizon T-L + , part 
of which has come from the past event horizon TiT and 
the rest from T~ . 

We remark that our "out" and "down" modes are not 
the same as those considered, for example, in 4J]. In 
[Zij , "out" and "down" modes are constructed from "in" 
and "up" modes by writing them in terms of Kruskal 
co-ordinates, taking their complex conjugates, and re- 
versing the signs of the Kruskal co-ordinates. This pro- 
cedure yields mode functions which are non-vanishing 
only on the left-hand-diamond of the Kruskal diagram 
(denoted region IV in Fig. [1}. However, the "out" and 
"down" modes that we have constructed (|2.35H2~3l)]) are 
non- vanishing on the right-hand-diamond of the Kruskal 
diagram (denoted region I in Fig.[T]). 



III. QUANTUM FIELD THEORY OF 
FERMIONS ON KERR 



Before we study in detail the definition of quantum 
states for fermions on Kerr black holes, we review the es- 
sential features of fermion quantum field theory in curved 
space, particularly stressing how this differs from the 
quantum field theory of bosonic fields. 

The first step is to select a basis of solutions of the 
Dirac equation (|2.9[) which are orthonormal with respect 
to the inner product defined in (|2.21[) and expand the 
fermion field in terms of this basis. Before promoting the 
coefficients in this expansion to operators, it is necessary 
to divide the mode solutions of the field equation into two 
sets: the expansion coefficients of one set will correspond 



to particle annihilation operators, and the expansion co- 
efficients of the other set will correspond to particle cre- 
ation operators. We will denote the modes in the first set 
as ip A and the ones in the second set as tp A . This division 
of the modes is not completely arbitrary: it must be the 
case that the particle annihilation and creation operators 
satisfy the canonical quantization relations. One usu- 
ally chooses the modes i/j a as being 'positive frequency 
modes' with respect to a chosen time-like co-ordinate r 
(that is, when they are Fourier-decomposed with respect 
to t they only contain positive frequency components) 
and the modes ip A as being 'negative frequency modes' 
with respect to the co-ordinate r. If the space-time has a 
globally time-like Killing vector d/dr, then the choice of 
'positive frequency' using the co-ordinate r is the most 
natural and corresponds to positive frequency modes also 
having positive energy. 

For a quantum scalar field (see App. [A] for Kerr space- 
time and [45| for rotating Minkowski space-time) , particle 
annihilation operators dA and creation operators a A are 
conventionally taken to satisfy the commutation relations 



a a, a a' 



3AA' 



(3.1) 



where [d, b] = ab — ba. The canonical commutation rela- 
tions (|3.ip have the consequence that the modes which 
arc the coefficients of the operators a a (that is, the 
equivalent of the modes ipt for fermions) have positive 
Klein-Gordon norm and the modes which are the coef- 
ficients of a A have negative norm (essentially, if it were 
the other way around, the sign on the right-hand-side of 
(|3.1[) would change, leading to an interpretation of a as a 
particle creation operator and d^ as a particle annihila- 
tion operator). For quantum scalar fields, the sign of the 
norm of the mode in general depends on the frequency, 
which restricts the possible choices of positive and nega- 
tive frequency modes as, respectively, coefficients of the 
annihilation and creation operators [45J . 

For a quantum fermion field, the particle annihilation 
and creation operators are conventionally taken to satisfy 
the canonical anti-commutation relations 



{d A ,d A ,} 



= 6 



AA' 



(3.2) 



where {d, b} = ab + ba. All the fermion modes defined 
in Sec. Ill CI have positive norm, independent of the fre- 
quency of the mode (the same is true for fermions in 
rotating Minkowski space-time [3^,[1^]). In other words, 
for fermion fields, both positive and negative frequency 
modes have positive norm. This means that, unlike the 
scalar case, positivity of the norm does not restrict the 
choice of positive and negative frequency modes as coeffi- 
cients of the annihilation and creation operators. There- 
fore we have rather more freedom in the fermion case to 
choose the modes tp A according to physical criteria, for 
example requiring the energy of a mode as seen by a par- 
ticular observer in a particular region of the space-time 
to be positive. 
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We now outline the construction of our fcrmion quan- 
tum field theory. We consider an orthonormal basis of 
modes of the form (|2.22j) . Next, we make an appropri- 
ate choice for the modes ip^ and ip^ (see the rest of this 
section for the physically relevant choices). The ip\~ and 
ip A then form a complete orthonormal basis of solutions 
of the Dirac equation. Note that the spinor modes ip A 
have the form (|2.22[) and are not the complex conjugates 
of the spinor modes ip A because of the complex function 

We therefore expand our fermion field ^ in terms of 
these basis spinors: 



1. 'Past' Boulware state \B } 

On X~ , it is natural to define positive frequency with 
respect to the Boyer-Lindquist time co-ordinate t, since 
this is the proper time for an observer at rest far from 
the black hole. A suitable set of modes having positive 
frequency with respect to t on T~ is 



where lu > 0, 



1 



(3.6) 



(3.3) 



Va 



2B i l(r) 2 SA(9) 



(3.7) 



where the sum is over the appropriate values of the quan- 
tum numbers A. We note that the expansion coefficient 
b\ is, unlike the case for a scalar field, not equal to a\. 
Quantization proceeds by promoting the field $ and ex- 
pansion coefficients cla and 6a to operators, with the lat- 
ter satisfying the anti-commutation relations 

{a A ,a A ,} = <W = {&A,&l}. ( 3 - 4 ) 

We interpret the operator a A as an annihilation operator 
for fermions, the operator a A as a creation operator for 
fermions, and 6a, b A as annihilation and creation oper- 
ators for anti-fermions, respectively. We note that the 
annihilation operator for a fermion is not the same as 
the creation operator for an anti-fermion. The vacuum 
|0) is then defined as that state which is empty of both 
fermions and anti-fermions: 



a A \0) =0 = 6 A |0). 



(3.5) 



and the "in" radial functions are given by (|2.32[) for L — 
+1, and by (j232"j) with xi? A O 2 Ra for L = -1. 

The 'past-Boulware' state \B~) 0, [H| is defined by 
expanding the quantum fermion field ^ in terms of the 
above "in" modes (|3.6[) plus a set of "up" modes with 
positive frequency with respect to t on the past event 
horizon H~ . From the form of the radial functions (|2.33[) 
near the past event horizon, the relevant frequency near 
the event horizon is not lo, but lu = lu — mfln instead. 
This is because the "up" modes should be written in the 
form 



i 



(3.8) 



where <p = ip — flnt is the azimuthal co-ordinate which 
co-rotates with the event horizon, and 



up _ ( iR A P (r)iSA(e) 
- { 2R A P (rhSA(e) 



(3.9) 



It is clear from the above construction that, as with scalar 
fields, the definition of the vacuum |0) depends crucially 
on the choice of the modes which are the coefficients of 
the operators cla- 

What is different about the fcrmion field, however, is 
that there is much more freedom in making this choice, 
which will be of fundamental importance for the rest of 
this section. 



A. 'Past' and 'future' quantum states 

The surface %~ Ul~ is a Cauchy surface for the right- 
hand-quadrant of the Kruskal diagram for Kerr (denoted 
by region I in Fig. [TJ, in other words for the space-time 
exterior to the event horizon. We begin by reviewing 
the construction of quantum states defined in terms of 
properties on this Cauchy surface, since this is uncon- 
troversial and can be performed for bosonic as well as 
fermionic fields. All the states we consider in this section 
are not invariant under simultaneous t — tp reversal. 



the "up" radial functions being given by (|2.33l) for L = 
+1 and by (|2~33)l with ii? A ^ 2-Ra for L = -1. For the 
modes in (13.81). we have 



di 



V7 = 



- ~ /UP 



(3.10) 



so that the natural choice of positive frequency for the 
"up" modes near H~ is uj > 0, reflecting the fact that 
an observer near the event horizon cannot remain at rest 
relative to infinity. 

The modes (|3.6p and (|3.8p form a complete orthonor- 
mal set and therefore, splitting the field into modes ip A 
and tp A and following the procedure outlined at the start 
of this section, we expand the quantum fermion field as 



(lid 



^7 + &r 



(3-11) 
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where we remind the reader that —A = {— uj,£, — m}. 
The expansion coefficients have become operators satis- 
fying the usual anti-commutation relations 

{a^ n/up , a£ /up } = 5 AA , = {^ l/up , S£ /up } , (3.12) 

where A' = {uj',i', m'} and <5aa' = 5 (uj — uj') Si^>S mim '. 
The 'past-Boulware' vacuum \B~) is then defined as that 
state annihilated by the a and b operators: 

Sfc\B-) = i%\B-) = 0, uj>0, 
a u A p \B~) =b u ^\B-) = 0, w>0. (3.13) 
I 



This definition of the 'past-Boulware' state is the same 
as for the bosonic case (modulo the subtleties in defin- 
ing the "u p" modes for bosons), and is the state con- 
sidered in [15|. It corresponds to an absence of parti- 
cles either coming in from X~ or emanating from the 
past event horizon TiT . However, this state is not a vac- 
uum state as seen at I + : it contains an outgoing flux 
of particles in the "up" modes where ujuj < 0, which is 
the Unruh-Starobinskh radiation [l5|,[l||. This 'quantum 
super-radiance' occurs even though fermions do not dis- 
play classical super-radiance (see remarks below (|2,38[) h 



2. 'Past' Unruh state \U~) 



Next we turn to the definition of the 'past-Unruh' state \U~) 0, [28|. The "in" modes (|3.6p are again chosen to 
have positive frequency with respect to Boycr-Lindquist time near I - . However, we now require the "up" modes 
(|3.8p to have positive frequency with respect to the Kruskal retarded time (that is, the affine parameter along the 
null generators of the past horizon [13]) near the past event horizon %~ . Using the Lemma in Appendix H of [HI, it 
can be shown that a suitable set of positive frequency modes is given by the following, for all values of uj [2l| : 



2 cosh 



\2T H J 



exp 



4Th 



^down, 



(3.14) 



where Th is the Hawking temperature of the black hole ()2.4|) . There is a subtlety in the definition of the ip A wn 
modes: these are obtained by taking the complex conjugate of the "up" modes and changing the sign of the Kruskal 
co-ordinates. The ^ own modes arc therefore not the same as our "down" modes f/;^ own formed from the radial 
functions (|2.36l) : the latter are non- vanishing on the right-hand-quadrant of the Kruskal diagram for Kerr (region I in 
Fig- HI): while the former are vanishing on the right-hand-quadrant of the Kruskal diagram and so do not need to be 
considered in detail. Similarly, a suitable set of modes having negative frequency with respect to Kruskal time near 
TiT is found to be, again for all values of uj: 



2 cosh 



2T H J 



exp 



■iky- 



exp 



4Th 



^down, 



(3.15) 



Further details of this construction can be found in 0, |44| . We therefore expand the quantum fermion field in terms 
of these positive and negative frequency modes as follows, where we work on the right-hand-quadrant of the Kruskal 
diagram only: 



E E 

£=1 m= — t 



duj 



duj 



2 cosh 



*P7 



exp 



4T ff 



(3.16) 



The 'past-Unruh' state \U ) is then defined as that 
state which is annihilated by the c and al operators: 

£%\U-) = a*£\U-) = 0, uj>0, 
c u A p \U-) = d A p \U~) = 0, all uj. (3.17) 

As with the 'past-Boulware' state \B~), the derivation 
above mirrors that for the bosonic fields (see, for exam- 



ple, Appendix B of 0), except that for fermions there 
are no difficulties in defining the "up" modes. The 'past- 
Unruh' state \U~) corresponds to an absence of particles 
incoming from X - , but, as we shall see in Sec. IIV1 the 
"up" modes from H~ are thermally populated. 

The 'past-Boulware' and 'past-Unruh' states defined 
in Sees. MI A II and MI A 21 are uncontroversial and well- 
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defined for quantum fields of all spins. Various expecta- 
tion values in these states have been computed for both 
fermionic and bosonic fields, see [H, [13, 0, HH, [H, |28[ 

Hi, Si, Si. 



3. CCH-state \CCH~) 

There is one further 'past' quantum state which can 
be defined. For the 'past-Unruh' state above, there is an 
absence of "in" mode particles but the "up" modes are 
thermalized with a thermal factor containing their natu- 
ral mode energy Co. One can define a further state, the 
Candelas-Chrzanowski-Howard (CCH) state [23], which 
we denote \CCH~) (see App. lA~2]l . In the \CCH~) 
state the "in" modes are thermalized as well as the "up" 
modes, using the natural mode energy w in the ther- 
mal factor for the "in" modes. In common with the 
other 'past' quantum states considered in this section, 
the CCH-state \CCH~) is not invariant under simulta- 
neous t — ip reversal. For bosonic fields, expectation val- 
ues in this state have been found to have good regularity 
properties (29|. 



4- 'Future' quantum states 

Following [28j], we could use "out" and "down" modes, 
defined from the radial functions (|2.35H2~3l))) and con- 
sidered in more detail in the next section, to define a 
'future-Boulware' state \B + ) which would correspond to 
an absence of particles from I + and H + . We do not con- 
sider this further in this article; instead, in Sec. IIII B"l we 
will define a state which is empty at both 1~ and I + . 

It would also be possible to define a 'future-Unruh' 
state \U + ) (28| by considering "out" modes with positive 
frequency with respect to time t at I + and "down" modes 
with positive frequency with respect to Kruskal time near 
H + . This state would have no outgoing particles at X + 
but the "down" modes would be thermally populated. 

In analogy with the 'future-Boulware' and 'future- 
Unruh' states above, we could also define a state \CCH + ) 
by thermalizing the "out" and "down" modes with their 
natural energies appearing in the thermal factors. We do 
not consider such 'future' states further in this paper. 



and W + . The existence of one of these two states is inti- 
mately linked with the existence of the other. 



B. A candidate 'Boulware' state 

For scalars and electromagnetic radiation, it is shown, 
respectively, in [28| and [23 that a 'Boulware' state, 
empty at both 1~ and I + cannot be defined (see also 
App. IA 2[) . Instead one has to consider the 'past- 
Boulware' \B~) (see Sec. IIII A ip and 'future-Boulware' 
\B + ) (see Sec. IIII A~4f states constructed in the previous 
subsection. 

However, we now show that for fermions the situa- 
tion is different. We have already defined a set of "in" 
modes Q3.6P which have positive frequency with respect 
to Boyer-Lindquist time t at 1~ . Similarly, a set of "out" 
modes, having positive frequency with respect to t at I + 
can be defined as follows: 



where u > 0, 



1 



out 

-nut imip I 'I A 



out _ ( iRT (OiSaW 

Va - { 2Rr t (rhS A (6) 



(3.18) 



(3.19) 



and the "out" radial functions arc given by (|2.35l) for 
L = +1 and by ((2351) with x i? A <-> 2 R\ for L = -1. 
Expanding the fermion field in terms of the "in" and 
"out" modes gives 



V-AeT + ^-a/a 



in f'm\ 



duj 



out ?outt 
A J A 



+ Ca/ 



(3.20) 



As with scalar fields (sec App. I A lj) . the "in" and "out" 
modes, while complete, are not orthogonal and so they 
do not give a suitable quantization basis. Instead, we 
write the "out" modes in terms of the orthogonal "in" 
and "up" modes, using the relations (|2.39|) : 



0a 



out 
A 



A OUt 



A ' 



down ^downj up 

- A A ^A " 



ndown 
B A 



A ■ 



(3.21) 



We now turn to the more subtle task of defining 'Boul- 
ware ' \B) [3 and 'Hartlc-Hawking' \H) @ states for 
fermions on Kerr. By a 'Boulware' state, we mean a 
state which is empty at both T~ and X + . By a 'Hartlc- 
Hawking' state, we mean a state which represents a ther- 
mal bath of radiation at the Hawking t emp erature of 
the black hole. It would be anticipated [25[ that such 
a 'Hartlc-Hawking' state, if it exists, would respect the 
symmetries of the space-time and be regular on both TiT 



noting that this transformation only involves ipA^ UP ano - 
not their complex conjugates (in contrast with the scalar 
case in the super-radiant regime, Eq. (|A8[) ). and is valid 
for all signs of w and oj. We define the modes ^ own 
similarly to -0a P in Eq. Q3.8P but using the radial functions 
i,2-R A own (given by (l2~36l) for L = +1 and by (l2~36)l with 
\Ra O 2-Ra for L = —1) instead of i,2-Ra P - The relations 
()3.2ip enable us to rewrite the expansion (|3.20p , now for 
the quantum fermion field, in terms of "in" and "up" 
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modes: 



oo I 



^y+^a/a 



inf 



d'jj 



(3.22) 



where 



-out /tout 



A' 

?inf ?inf . ibutf ,iout 

/a — /a T Ja ^A j 



3 up 



-out^out 



A J 



/A Pt = /A Utt SX ut - (3.23) 



Again, the expansion coefficients have become operators 
satisfying the usual anti-commutation relations, and then 
we define our candidate 'Boulware' vacuum \B) as that 
state annihilated by the e and / operators: 



\B) 



Ia\B) 



\B) 



fr\B)=o, ^>o. 



(3.24) 

Of course, the fact that we have defined a candidate 
'Boulware' state does not mean that this state is regular 
or Hadamard (anywhere), or, indeed, physically relevant. 
However, it is worth stressing that, in the fermion case, 
we have been able to progress rather further with the 
definition of a candidate 'Boulware' state than is possible 
with bosonic fields. Note that at this stage we are not 
making any claims whatsoever as to the regularity of the 
state \B); instead we are simply commenting that, at this 
stage, our definition seems reasonable. In Sec. IIV CI we 
will compute some differences in expectation values for 
observables between two states, including the state \B), 
which will provide concrete evidence for the existence of 



this state and its regularity, at least on part of the space- 
time exterior to the event horizon. 



C. A candidate 'Hartle-Hawking' state 

The Kay-Wald theorem [25| proves that in essentially 
any globally-hyperbolic and analytic space-time with a 
bifurcate Killing horizon there can exist at most one 
Hadamard state which is regular everywhere and respects 
the symmetries of the space-time. The theorem further 
proves that, if such a state exists, then it must be a ther- 
mal state. Importantly, Kay and Wald show that such 
a state does not exist for scalar fields on Kerr space- 
time. Therefore there cannot exist a 'Hartle-Hawking' 
state which is regular everywhere outside the event hori- 
zon and on both T~L~ and H + . While the Kay-Wald re- 
sult is proved formally only for scalar fields, one could 
anticipate that it is valid for fields of higher spin, in- 
cluding fermions. Of course, the Kay-Wald result is a 
non-existence theorem, and it may be possible, for ex- 
ample, to have a state which respects the symmetries 
of the space-time but is not regular everywhere. For 
scalars, Frolov and Thorne Q have used the ^-formalism 
to construct the so-called FT-state (see Add. I A 2"]) . which 
respects the symmetries of the space-time but unfortu- 
nately is well-defined only on the axis of rotation [28| . In 
this section we will construct a state which possesses the 
symmetries of the space-time, before undertaking some 
numerical computations in Sec. IIV CI to investigate its 
regularity properties. We emphasize that we do not need 
to use an analogue of the ^-formalism for defining this 
state for fermions. 



1. 'Hartle-Hawking' state \H) 



To define a state which has the potential to be regular on both W~ and H + , we seek modes which have positive 
frequency with respect to Kruskal time near both H~ and H + . In Scc. lIII Al we have already constructed positive and 
negative frequency modes with respect to Kruskal time near rl~ (see p.!4[) and (|3.15|) respectively). By a similar 
method, a suitable set of modes having positive frequency with respect to Kruskal time near H. + is found to be, for 
all Co: 



2 cosh 



2T H ) 



exp 



f— 



^,dow„ 



exp 



(3.25) 



and a suitable set of modes having negative frequency with respect to Kruskal time near Ji + is found to be, for all a): 

(3.26) 



2 cosh ( ) 
\2ThJ 



exp 



^down 



CXP fe)^ 



In (|3. 25113. 26]) . as in p. 14113. 15]) . the rp^ modes are defined by taking the complex conjugate of the "down" modes 
and changing the sign of the Kruskal co-ordinates. The modes are therefore not the same as our "up" modes 
ip^ , and vanish on the right- hand-quadrant of the Kruskal diagram (region I in Fig. [!}. As in Sec. 1111 A 21 we do not 
need to consider them further. 
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We therefore expand our fcrmion field on the right-hand-quadrant of the Kruskal diagram in terms of the modes 
(j3~T3H3"J5l |3~23H3~2S)) to obtain: 



* = EE 



2 cosh 



2T H J 



cxp 



9 a + ex P 



4Tff / 



\ rupt 

A 



^down 



exp 



AT, 



H 



3A° wn + cxp 



AT H 



^downf 



(3.27) 



The "up" and "down" modes are not orthogonal so do not form a good quantization basis. We use the relations (j2.39[) 
to write the "down" modes in terms of "in" and "up" modes (we could equally well write the "up" modes in terms 
of "out" and "down"), obtaining, for the quantum fermion field: 



oo I 



£=- rr 

— 2 



exp 



i 

(— 

\AT H 



2 cosh 



V2tJ 



3a + exp 



i>7 



>t 



exp 



( JL \ ~ u p 



9 a + ex P 



AT H 



A 



AT H 



(3.28) 



where 



~up ~ U P i -down /(down 

9 A — 9 A + 9 A A A i 

h7 ] = + h^AjT 



-in -down ndown 
9 a = 9 A B A 



^, in t 



(3.29) 



We then define our candidate 'Hartle-Hawking' state \H) 
as that state which is annihilated by the g and h opera- 
tors: 



9 [ a\H) 



h%\H) 



97\H) 



h u A p \H)=0, VCj. (3.30) 



As with our candidate 'Boulware' state in Sec. IIII Bl we 
cannot at this stage make any claims as to the regular- 
ity or properties of our candidate 'Hartle-Hawking' state. 
However, we are encouraged by the fact that we have 
been able to proceed this far for fermions (the corre- 
sponding construction for bosons fails due to the super- 
radiant modes and the need to use positive norm modes) . 

Further evidence that our candidate 'Hartle-Hawking' 
state | H) may be regular on at least part of the space- 
time exterior to the event horizon is provided by consid- 
ering the simpler situation of a rigidly rotating thermal 
bath in flat space, as, at least close to the event hori- 
zon, it is expected that a 'Hartlc-Hawking'-like state on 
Kerr space-time should represent a thermal bath of ra- 
diation rotating rigidly with the angular velocity of the 
event horizon fi^-. For a rigidly rotating thermal bath 
of scalar particles in flat space (50l |. the quantum state 
is ill-defined everywhere. The situation for fermions in 
flat space is rather different [46j]. It is possible to define 
a state in flat space which is regular inside the specd- 
of-light surface Sl but diverges on Sl (and, presumably, 
outside Sl as well) . These results indicate to us that our 
'Hartle-Hawking' state on Kerr should be defined and 
regular, at least sufficiently close to the event horizon. 

One final comment is in order in this section, namely 
can we construct an analogue of the Frolov-Thorne state 



for fermions? We will see in Sec. IIVI that expecta- 
tion values of operators in the Frolov-Thorne state for 
fermions can easily be defined and turn out to be iden- 
tical to those for our candidate 'Hartle-Hawking' state 
\H). We will therefore conclude that our new state \H) 
is indeed the fermionic analogue of the Frolov-Thorne 
state. 



2. An alternative vacuum state \B) 

For further comparison with both the scalar and 
fermion field results for a rigidly rotating thermal bath 
in flat space-time 0, it is helpful to have, for 
Kerr space-time, an analogue of a 'vacuum' state which 
is defined within the speed-of-light surface (our candi- 
date 'Boulware' state for Kerr space-time, constructed in 
Sec. IIIIB( is not helpful in this regard because it is de- 
fined with respect to infinity and we suspect that it may 
not be regular all the way down to the event horizon). 
To do this, we expand the fermion field in terms of "up" 
and "down" modes with Cj > 0: 



duj 



du 



^/down-down i jdown -downf 
Va X A + V-A Va 



(3.31) 



As with the candidate 'Hartle-Hawking' state \H) (see 
Sec. IIII C"Tj) . we write the "down" modes in terms of the 
"in" and "up" modes, and then, promoting the resulting 
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expansion coefficients to operators, we find 



* = E E 



P— 1 ra- 
ti, Q 



dtu 



^7 + tfw* 



(3.32) 



We then define yet another vacuum \B) as that state 
annihilated by the x and y operators: 

xl\B) = ff£\B) = xf\B) = y u A P \B) = 0. (3.33) 



number current operator J M and stress-energy tensor 
operator in each of the states defined in Sec. IIII1 
namely 'past-Boulwarc' \B~) (Sec. MI A lj) . 'past-Unruh' 
\U~) rSec. lIII A2|) . the CCH-state \ CCH~ ) (Sec. IIII A 3ft . 
our candidate 'Boulware' \B) (Sec. IIII Bj) . our candidate 
'Har tle-Haw king' \H) (Sec. IIII C 1)) and the state \B) 
( Sec. IIII C 2j) . Unfortunately, rcnormalization of all these 
quantities on Kerr space-time remains an intractable 
problem, and therefore our analysis is limited to find- 
ing the differences in expectation values between two of 
the above states. 



Once again, at this stage we make no claims as to the 
regularity of the state \B), merely that the definition 
above seems reasonable. In particular, we should em- 
phasize that the state \B) is not a candidate for a state 
on Kerr analogous to any of the standard Schwarzschild 
black hole states (Boulware, Unruh or Hartle-Hawking). 
We have introduced this state solely to aid the interpre- 
tation of the state \H) in Sec.[V] We expect that the state 
\B) will approximate a rigidly rotating vacuum state with 
the same angular speed as the event horizon, analogous 
to the fermionic rotating vacuum in flat space [46j . 



IV. EXPECTATION VALUES OF 
OBSERVABLES 

We now turn to the computation of the expectation 
values of various observables in the quantum states de- 
fined in Sec. IIII| in particular to investigate the prop- 
erties of our candidate 'Boulware' and 'Hartle-Hawking' 
states. We are interested in expectation values of the 

I 



A. Observables 

The simplest non-trivial fermion operator to study is 
the number current J M , given as a quantum operator by 



■>' = \ 



(4.1) 



In (|4.1[) , the commutator is understood to act only on 
the operators in \P and not on the spinor mode functions, 
which keep the order ipj^ijj so that expectation values of 
do not have any spinor indices. Physically, expecta- 
tion values of the operator J 1 , i = 1,2,3 count the flux 
of particles (that is, flux of fermions minus flux of anti- 
fcrmions) in a particular direction and the expectation 
value of J* counts the particle number density (again of 
fermions minus that of anti- fermions). Note that these 
quantities will not be zero in general because the black 
hole emits fermions preferentially in the southern hemi- 
sphere and anti-fermions in the northern hemisphere [l8l - 



The expectation values of J M (|4. 1[) in each of our states of interest can be written in terms of the classical number 
current f|2 . 20[) acting on individual modes as follows: 



(B~\J^\B- 



(CCH-\P\CCH- 



(B\J»\B) 



1 oc £ 

\y. e 

f=i m=—£ 



(u-\.hu-) = -J2 E 



1=1. m= 



E E 



(b\p\b)=- y: e 



oo e 



t=\ m=-£ 
oo I 

EE 



{H\P\H) = -Y, E / ^ tanh 



duj a' 



>Ja 



duj tanh 



duj 



duj 



du> tanh ( | j^ p,/i 



\2ThJ 



2T h J Ja 



J A + 3 A 



2T H ) 



■m.LL , up, a 



du tanh ) 



3 A + J A 



(4.2) 
(4.3) 
(4.4) 
(4.5) 
(4.6) 
(4.7) 
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where 



• in/up, u -7-in/up , in/up -rin/ u P , in/up 



(4.8) 



We can also write down the expectation values of J M for the analogue of the FT-state \FT) Q: 



(FT\J»\FT) = lt± {£ *> tanh ( J-) ft" + [ to tanh ( J-) 



(4.9) 



noting that this differs from (|4.4j) in the thermal factor for the "in" modes. At first sight, it looks like (|4.9|) differs from 
the expecation value for our candidate 'Hartle-Hawking' state \H) (|4.6[) in the integral over the "in" modes. However, 
it can be shown that the expectation values of the fermion current (and stress-energy tensor) for the FT-state \FT) 
and our candidate 'Hartle-Hawking' state \H) are in fact equivalent, so that, for all practical purposes, the fermion 
FT-state is the same as our state \H). 



Writing out the spinor mode functions explicitly in terms of the radial and angular functions using p. 61 13. 8|) . the 
classical mode contributions to the components of the current J A = V'aT^'V'A are (where we omit the "in/up" mode 
labels as these formulae apply equally well to all modes): 



7* 
J A 



T r 
■J \ 



J" 



TV 

J A 



1 



47r 2 A£sin 
L 



{zaix/Asin^ [xi? A2 i?A - iRa 2 R* a ] iS A2 S A - {r 2 + a 2 ) [UR A \ 2 x 5 a + | 2 i? A | 2 2 S 2 A ] } , (4.10) 



47r 2 Esin( 
L 



I i -Ra 1 2 iS\ — |2-Ra| 2 2 S\ 



47r 2 VASsin 
1 

~47r 2 AS]sin 2 6> 



iR A 2Ra + 1RA2RX} iSa2S\, 
a {^Lv / A hRhRA - 1 R A2 R* A ] iS A2 S A -asinO [\iR A \ 2 1 S 2 A + \ 2 R A \ 2 2 S 2 A ] } 



(4.11) 
(4.12) 
(4.13) 



From the above expressions, using the symmetries (|2.28H2"T2"9"|) . we find 



A 



3 a 



Ja 



J A = 



47r 2 ASsin6» 
L 

47r 2 Ssin 9 
L 



iRa\ 2 -\2Ra\ 2 



[iS a - 2 S A ] 



\iRa\ 2 + \2Ra\ 2 



\iS a 



iSa] ; 



7r 2 a/AS sine 
a 

47r 2 ASsin( 



; — -3? (li?A 2-Ra) iSa 2 S A , 



l-Ra| 



iRa\ 



A - 2S A ] , 



(4.14) 
(4.15) 
(4.16) 
(4.17) 



where again we have omitted the superscript in/up because the above expressions apply equally well to "in" and "up" 
modes. The expressions (|4.10H4~T7j) depend explicitly on L. In view of our comments in Sec. Ill CI regarding how to 
obtain the expressions for L = —1, we note that if one uses the boundary conditions as written out in Eqs. (|2.32| - |2~3"3"|) . 
then Eqs. (|4.14fl4T7j) are already valid directly for L = +1. On the other hand, for L = — 1, if one chooses to continue 
using the boundary conditions Eqs. (|2.32||2~3"3)) . then Eqs. (|4.14ffl~TT)) arc valid by setting L = — 1 and also swapping 
iR A 2Ra in these latter equations. 

In the absence of a framework in which to perform computations of renormalized expectation values on Kerr space- 
time, in this article we study differences in expectation values in two different states. The particular differences on 
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which we focus are: 

oo I 

(>) U ~~ B ~ = (U-\P\U-) (B-\J»\B-) = - E E 
(P)Cch--b- = (CCH-\J»\CCH-) - (B~\P\B 

oo I 

=-E E 



die 



l + e T H 



3 A ' 









/•oo 








1 + e T « 




/ da) 


1 + e T « 








Jo 





(fl|>|B)-<B-|>|B-) = ig E 



oo I 

-E E 



(B-\J»\B- 



(J") 



f=i m= - 

(H\J"\H) - (B\Ji*\B 









POO 






{f ** 


1 + e T « 




1 duj 


1 + e T » 










Jo 







E E 



1 



1 r . 



ti\H-B 



{H\J^\H) - (B\J^\B) 

oo I 

-E E 

£=i TO=- 



/•OO 




-1 r 


/ do) 


1 + e T « 
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in,/i i -up,/i 
J A + .7 A 



.7 A + J A 



E E 



do; 



m,/i , -up, fi 



Ja 



3a 



(4.18) 

(4.19) 
(4.20) 

(4.21) 

(4.22) 

(4.23) 



in terms of which all other differences in expectation values can be computed. In (|4.18hl4~2"3)) , we have introduced the 
notation (0) A ~ B = (A\0\A) — (B\(D\B) for the difference in expectation values of the operator O in the states \A) 
and \B), and we shall use this notation for the remainder of the paper. 

I 



The main observable of interest is the expectation -0a): 
value of the stress-energy tensor operator . As a quan- 
tum operator, T^ u is given by aT!^/ up 



'■yin/up in/up . "jin/up , in/up 



} , (4-24) 



—in/up 
V^A 

— in/up 



7^i n/UP 



7/x^, 



i/up 



(4.26) 



where, as with the number current operator, the commu- 
tators are understood to act on the operators in # and 
not on the spinor mode functions, which retain the order from those" 
i/j"f f j,tjj. The expectation values of in our states of in- 
terest take the form (I4.2H4.7I). but with the mode contri- 



The expressions for aT™J up and a^J?/ up are extremely 
lengthy, so we relegate them to App. [C] As with the 
number current we are interested in differences in expec- 
tation values of T^ v between two states; the key ones are 
in (|4.18fl4~23)) , and all other differences can be computed 



butions to the current j)^ up ' M replaced by the quantities 
a£™/ up , where 



B. Numerical method 



At 



in/up 



/ ' ' ' 



= - A t 1 



in/up 



/ ' ' ' 



in/up 



Here we address the challenge of computing the dif- 
ferences in expectation values of quantum states nu- 
(4.25) merically, by evaluating their mode sum representations 
(14.181 l4~23l) . Computing a typical example (X) is not 
a trivial task, for a number of reasons. Firstly, (X) is 
and aT^J uv is the classical mode contribution to the a function of radial and angular coordinates r, 9, and so 
stress-energy tensor (that is, (|2.19|) with replaced by must be evaluated on a representative grid of points. Sec- 
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ondly, for each grid point, (X) is computed from a double 
sum and an integral over frequency, 



(*>M) = E E 



X im (Lo;r,9)dto. (4.27) 



Thirdly, the integrand X£ m (to; r, 9) is computed from ra- 
dial and angular functions, i,2-RA^ UP ( r ) an< ^ i.2Sa(0), 
which are obtained from the numerical solutions of or- 
dinary differential equations (|2.25[ 12 . 2T[) , with appropri- 
ate boundary conditions (|2.321 |2~3"3"|) . Finally, (X) is not 
necessarily finite and well-defined in some subregions of 
the (r, 9) plane, for example: at the horizon, inside the 
stationary limit surface, or outside the speed-of-light sur- 
face. 

We first outline our method for finding the summands 
Xi m (u)\r,6) (that is, computing the radial and angular 
mode functions), before turning to the computations of 
the mode sums and related convergence issues. 

1. Mode functions 



To compute the angular eigenvalues A (see (|2.27p and 
cigenfunctions i,2Sa((9) we applied the spectral decompo- 
sition method described in j53[, in which i^Sa^) is ex- 
pressed as a series of spherical spin-half harmonics. This 
approach leads to a three-term recurrence relation for the 
coefficients of the series, and the convergent solution may 
be found via the method of continued fractions (see, for 
example, (lij). We checked our results by implementing 
an alternative three-term relation given in [55|. Typi- 
cal angular functions are shown in Fig. [3] (a) , for I = | , 
m = i and a range of values of aio. The plot shows that 
the symmetries (|2.291 12.30[) are satisfied by our numerical 
angular functions. 

The "in" and "up" radial functions i t 2R^ up (r) are 
found from numerical solutions of (|2.25[ 12.270 subject 
to boundary conditions (|2.32[ I2.33[) . To compute these 
modes we made use of generalized series expansions, in 
i — r/j at the horizon (for the "in" modes) , and in powers 
of 1/r at spatial infinity (for the "up" modes), as initial 
data for a Runge-Kutta integrator. The method closely 
follows the steps described in 2l|. Typical radial func- 
tions for the "in" and "up" modes are shown in Fig. [3] 
(b) and (c). 



w by a thermal factor (exp(w'/Tff ) + 1) _1 (where lo' G 
{to, Co}) which acts as a high-frequency cut-off (see Fig.rj|. 
Hence, for a given £, m (and < r < oo), the integral 
over frequency is finite. Furthermore, for a given £, the 
sum over m is finite. This leaves the infinite sum over £ 
as the only possible source of divergence. 

To perform the integral over frequency in ()4.27[) (for 
each r, 9, £, m) we first sampled the integrand over a uni- 
form grid of points across the domain of integration, af- 
ter replacing the infinite upper limit with a finite cut- 
off (if necessary), typically lo = max(0,mf2£r) + IOTjj + 
0.2M -1 . Then we interpolated the data with a cubic 
spline, resampled, and applied Simpson's rule to find the 
integral. The finite sum over m was straightforward to 
perform, whereas the infinite sum over £ required more 
consideration. We examined the contribution of the indi- 



vidual £- modes X/, and the truncated sum, with 



set 



to be a large value (typically ^ max ~ 20). The magnitude 
of these quantities gave an indication of convergence, as 
can be seen in Fig. U) 

In Fig. 5] we plot a typical integrand 



1 + exp ( — 



-i -l 



Ain) 



(4.28) 



(where the expression for tgg in terms of the radial and 
angular functions is given in (|C22p ) as a function of fre- 
quency lo, for "in" modes with ^<£<^,m>£ — 1, in 
the special case a = a ~ 0.910M. It can be seen in Fig.fJ] 
that modes with m = £ make the dominant contribution 
to the mode sum. For each fixed £, m, the integrand as 
a function of lo is strongly peaked at a particular value 
of to and the rapid convergence of the integral over lo can 
be seen. The location of the peak moves to higher values 
of lo as £ increases. In Fig. [4] (a), the magnitude of the 
peaks is decreasing very rapidly as £ increases past £ = | , 
indicating that the sum over £ is convergent in this case. 
In Fig. 2] (b) it is less clear whether the sum over £ is 
convergent or not, although the magnitude of the peaks 
of the integrand is decreasing at larger to. In Fig. @] (c) 
the peaks are still steadily increasing and the sum over £ 
does not appear to converge. 

A key part of our analysis is to determine whether 
or not the expectation values (X) are finite, so we con- 
clude from Fig. @] that a more sophisticated analysis of 
the mode sum convergence is required. If the terms 
in the sum are absolutely convergent, in the sense that 
lim^oo \X e /X e _i\ < 1, where 



2. Mode sums 

If (X) is not finite, then we would expect its mode sum 
representation to be divergent. To see how the divergence 
may arise, let us consider the ingredients in (|4.27|) . The 
mode functions -RaM and S\(0) are finite for r/j < r < 
oo. In cases where the frequency integral is taken over 
a semi-infinite domain (that is, when w max — > oo), the 
integrand in the frequency integral is suppressed at large 



X, 



E 



Xi m (Lo;r, 9)du>, 



(4.29) 



then the sum is clearly finite and well-defined. Con- 
versely, if the sum is not absolutely convergent then (X) 
may be ill-defined (at the very least, poorly represented 
by a sum over modes). Hence we may apply a simple ra- 
tio test to give an indicator of convergence, by examining 



Re = \Xi/Xi-i\ , 



(4.30) 
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(c) Radial functions : "up"-modes 



, a/M = 0.0 
, a/M = 0.5 
, a/M = 0.0 
, a/M = 0.5 



1.5 
1 

0.5 


-0.5 
-1 
-1.5 



r.\ '.. w \ 



Re 
Re 
Re 
Re 



lR|m, 
1 R lm: 
2 M lm. 



1 \ 
/I \ 



I \ 

I \ 

I \ 

I .'-I 



i \ ■ 
■. / i ,i 



/ \ 
/ \ 
/ ,\ 

I ■ V 



/ .'A 



\ •■\ t M ' 'i 
'••V ' 1 • V 

\ (•>■> : A /\ .• / 

V \../\ ^' K' 



'■■ i -. ^h?*- w^-if"- ■\\^ft n 



\ / 
\ / 



I = m = 3/2 
Mco = 0.4 
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r./M 



FIG. 3. Examples of typical angular and radial mode functions. Plot (a) shows the spin-half spheroidal harmonics iSa and 
2Sa for 1= |, m — i for a range of spheroidal couplings aui = —2,— 1.5, ...,1.5, 2. The symmetries (|2.29l 12.301) are apparent. 
Plot (b) shows the radial functions \Ra and 2Ra for the "in"-modes, defined by boundary conditions (|2.32|) . for I = m — 3/2, 
Muj — 0.4 and two cases: a/M = [dashed lines] and a/M — 0.5 [dotted lines]. Plot (c) shows the radial functions for the 
"up"-modes, defined by boundary conditions H2.33|) . with the same parameters. 
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FIG. 4. Frequency integrals and mode-sum convergence. These plots show a typical integrand, [1 + ex.p(G} /Th)\ 1 tffi 



(where 

t ee is defined in Eq. ljG22]l). as a function of frequency tv, for the "in" modes with \ <£<^-,m>l — 1, in the special case 
a = a ~ 0.910M. The integrand is evaluated on the equatorial plane (9 = n/2) at r = (a) 1.6M, (b) 1.8M and (c) 2.0M. The 
mode sum in cases (a) and (b) appears to be convergent, whereas the sum in the case (c) does not seem to converge. We note 
that the speed-of-light surface intersects the equatorial plane at r = 2M in this case, so plot (c) indicates that this particular 
mode sum will diverge outside the speed-of-light surface. The physical implications of this result are discussed in Sec. II V CI 
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as a function of r, 9. In Sec. IIV C II wc plot pp (for a 
large but finite value of I ~ 20) as a function of r, 9 to 
distinguish between divergent regions (where pp > 1) and 
convergent regions (where < 1). 



3. Validating our numerical results 

We validated our implementation with a few simple 
consistency checks. First, to test the radial functions, we 
numerically computed the Hawking flux using Eqs. (9)— 
(10) in [17|, and we verified that it matched the values 
given in Table I of [17] . Next, we considered an expres- 
sion for the energy flux as a function of angle, given by 
Eq. (2.12b) in [lgf , d 2 E/r 2 dfldt = A/E(f tr ). A subtlety 
here is that it is difficult to compute the flux for the Un- 
ruh state \ U~) directly (due to the lack of a large-w cutoff 
in the modal expressions (|4.3[l ). but rather easier to com- 
pute the flux for the state difference U~ — B, which may 
be found from mode sums (|4.18j) and (|4.20j) . The 'Boul- 
ware' state \B) is expected to be empty at infinity, and 
hence (asymptotically) the fluxes should be equivalent. 
In summary, wc computed 



lim 2ttA / (T tr 



U~-B 



sin 9 d9. 



(4.31) 



at large radii, and confirmed that it asymptotes to the 
correct energy flux as r — > oo, as given in Table I of firj ] . 



C. Numerical results 

We now present a selection of our numerical results, 
obtained using the methodology outlined in the previous 
subsection. First we examine where the quantum states 
defined in Sec. IIIII are regular, before turning to other 
physical properties of these states. 



1. Regularity of guantum states 

The first key question we wish to address is whether 
the quantum states defined in Sec. IIIII are regular outside 
the event horizon of a Kerr black hole. We begin, in 
Fig. [SJ by plotting the ratio pp (|4.30|) of successive terms 
in the f-sum for the differences in expectation values of 
the stress-energy tensor component Tgg given in (|4.18h 
I4.23[) . The component Tgg was chosen for this analysis 
because if the stress-energy tensor is regular in a freely 
falling frame crossing the event horizon (or stationary 
limit surface, or speed-of- light surface), then it must be 
the case that this component of the stress-energy tensor 
is regular (28j . 

Fig. [SJ shows the ratio p t (|4.30p . plotted for £ ~ 20, 
as a function of r, 9, with x = rsin# and z = rcos9. 
In Fig. [SJ the axis of rotation of the black hole is a 
vertical line through the centre of each diagram, and 



the equatorial plane a horizontal line through the centre 
of each diagram. The green dotted line is the speed- 
of-light surface; the purple dotted line the stationary 
limit surface (throughout this section we use the value 

a = ao = M-i/2 [\/2 — l] for which these two surfaces 

touch in the equatorial plane). The black circle denotes 
the region inside the event horizon. Divergent regions 
(where pp > 1) are blue and convergent regions (where 
pi < 1) are yellow. 

We consider first the uncontroversial 'past-Boulware' 
\B~) and 'past-Unruh' \U~) states, defined in 
Sees. MI A II and IIII A 21 respectively. From Fig. [SJ (a), 
it can be seen that the expectation value (Tgg) u ~ B is 
regular everywhere outside the event horizon, including 
inside the ergosphere and outside the speed-of-light sur- 
face. This is in agreement with numerical results for this 
expectation value for spin-1 fields [29j]. As will be dis- 
cussed in more detail in Sec. [V] we expect that both the 
\U~) and \B~) states will be regular everywhere outside 
the event horizon, and so, to examine the regularity of 
other states, it will be useful to consider the expectation 
values of those states relative to either \U~) or \B~). 

Next we turn to the state \CCH~) defined in 
Sec. IIII A 31 (27j . In Fig. [SJ (b), we can see (again in 
agreement with similar calculations for spin-1 fields [29[) 
that the expectation value (Tgg) CCH ~ B is regular ev- 
erywhere outside the event horizon, including inside the 
ergosphere and outside the speed-of-light surface. 

The next state to be considered is our candidate 'Boul- 
ware' state \B), defined in Sec. IIII Bl Fig. [SJ (c) shows 
that the expectation value (Tgg) B ~ B is regular every- 
where outside the stationary limit surface, but diverges 
inside the ergosphere. 

Finally, we consider our candidate 'Hartle-Hawking' 
state \H), defined in Sec. ImCl Firstly, in Fig. [SJ (d) 
we see that the expectation value (Tgg) H ~ B is regu- 
lar everywhere outside the event horizon and inside the 
speed-of-light surface (including the ergosphere), but di- 
verges on and outside the speed-of-light surface. Fig. [SJ 
(e) shows that the expectation value (Tgg) H ~ B diverges 
inside the ergosphere and outside the speed-of-light sur- 
face, but is regular between the stationary limit surface 
and speed-of-light surface. From Fig. [SJ(f), we see that 
the expectation value (Tgg) H ~ B also diverges outside the 
speed-of-light surface but is regular inside it, including 
inside the ergosphere. 

To further elucidate the behaviour of the states \H) 
and \B), in Fig. [SJ we plot the ratio pp, (|4.30[) for the 
expectation values 

(fgg) H - u ~ = (H\f ee \H) - (U-\f ee \U-), 

(4.32) 

{fee) § - B ~ = (B\fgg\B) - (B-\fgg\B-). 

(4.33) 

Comparison of Fig. [SJ (a) and Fig. [SJ (d) leads us to con- 
clude that the state \H) is regular between the event 
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FIG. 5. Ratio test to examine the divergence/regularity of expectation values of quantum states, for the stress-energy tensor 
component (Teg). The differences in expectation values for the six states defined in Eqs. (|4.18llT23l l are considered. In 

particular, these are: (a) {f eg } u ~- B ~ , (b) {f 0B ) CCH ~ - B ~ , (c) {T eg ) B - B ~ , (d) {T ee ) H - B ~ , (e) {f ge ) H - B , (f) <T w )*-f In 
each case, the ratio pi (|4.30|l is plotted for t ~ 20 as a function of r, 9, where z = rcos# and x — rsin#. The axis of rotation 
of the black hole is a vertical line through the centre of each diagram, and the equatorial plane a horizontal line through the 
centre of each diagram. The green dotted line is the speed-of-light surface; the purple dotted line the stationary limit surface 

(we use the value a = ao = M-\J 2 [s/2 — l] for which these two surfaces touch in the equatorial plane). The black circle is the 
region inside the event horizon. Divergent regions (where pi > 1) are blue and convergent regions (where pi < 1) are yellow. 
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(a): H - U" 



(b): B - B 





FIG. 6. Ratio test to examine the divergence/regularity of expectation values of quantum states, for the stress-energy tensor 
component (Tee)- The differences in expectation values for the states defined in Eqs. (14.321 [4733]) are considered. In particular, 
these are: (a) (Tee) H ~ U , (b) {Tee) B ~ B ■ The structure of the plots follows that in Fig. \5[ and the same parameters are used. 





Aj° 





FIG. 7. The expectation values (J^) 0011 ~ B for components of the fermion current, multiplied by A (|2.2[1 . The expectation 
values have been computed using (|4.14fj47T7)l with L = +1 (for L = — 1 the components have the same magnitude but the 
opposite sign). The expectation values are plotted on the vertical axis as functions of (r, 6), with x = rsin# and z = rcosd. 
In the horizontal plane, positive values are shaded in red, while blue denotes negative values. We use the value a = ao = 

Mi/2 [v2— l] for the rotation parameter of the Kerr black hole. 
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horizon and the specd-of-light surface, but divergent on 
and outside the speed-of-light surface. The divergence 
inside the ergosphere in Fig. [5] (c) is coming from the di- 
vergence of the state \B) inside the ergosphere, which can 
be seen in Fig. [5] (c) . From Fig. [6] (b) we conclude that 
the state \B), like the state \H), is regular between the 
event horizon and the speed-of-light surface but diverges 
on and outside the speed-of-light surface. 

Thus far we have restricted attention to the expecta- 
tion value of the component Tgg of the stress-energy ten- 
sor. While a divergence in this component is sufficient to 
render the whole stress-energy tensor divergent [28| , the 
regularity of (Toe) does not guarantee the regularity of all 
components of the expectation value of the stress-energy 
tensor, particularly at the event horizon. We there- 
fore consider the expectation values (J fM ) CCH ~ B (see 
Fig. and (T^) CCH ~~ B ~ (sec Fig. EJ) and for all com- 
ponents of the fermion current and stress-energy tensor 
(Fig. [5](b) implies that the component (f gg ) CCH ~ B is 
regular everywhere outside the event horizon). The ex- 
pectation values (T^ U ) CCH ~ B have also been studied 
in detail for quantum electromagnetic fields [29[ . 

From Fig. [7J the components of the expectation val- 
ues of the fermion current J M are all regular outside the 
event horizon (the regions shown in Figs. [7H5] include the 
ergosphere and part of the region outside the speed-of- 
light surface), and diverge on the horizon. Furthermore, 
all components apart from (J^) ' ' 11 ~ B flip sign under 
the mapping 9 — >• ir — 6 (which corresponds to z — > —z). 
This is due to the preferential emission of neutrinos in the 
southern hemis phe re and anti-neutrinos in the northern 
hemisphere [HH, U HI ■ 

From Fig. [5J all ten components of the stress tensor ex- 
pectation values are regular everywhere outside the event 
horizon, but diverge on the event horizon, with the ex- 
ception of the (t6) and (6(p) components, which appear 
to be regular on the horizon. These two components 
are much smaller than the others but are not identically 
zero. In [28[ it is shown that for scalar fields the (tO) 
and (0(p) components of the renormalized stress-energy 
tensor vanish due to the properties of the scalar mode 
functions; however this is not the case for gauge bosons 
[29| nor fermions, as seen here. From Fig. [8j it can be 
seen that all the components of the stress-energy ten- 
sor are symmetric under the mapping 9 — > it — 9 (which 
corresponds to z — >• —z) apart from the (t9), (rff) and 
(6 if) components, which flip sign under this mapping (as 
would be expected). 

Bringing together our results in this subsection, we 
conclude that the states \B~). \U~) and \CCH~) arc 
regular everywhere outside the event horizon. In anal- 
ogy with the situation for Schwarzschild black holes, we 
expect that the 'past-Boulware' state \B~) is divergent 
on both the future and past event horizons and that the 
'past-Unruh' state is regular on the future horizon % + 
but diverges on the past horizon 7-L~ . Accordingly, we 
conjecture that the state \CCH~) is regular on both the 



future and past event horizons. Of course, a full com- 
putation of the renormalized stress-energy tensor in this 
state would be necessary in order to verify our conjec- 
ture. Assuming these properties of the \B~) state, we 
deduce that the states \H) and \B) diverge on and out- 
side the speed-of-light surface but arc regular between 
the event horizon and the speed-of-light surface. Finally, 
we have evidence that the state \B) diverges in the ergo- 
sphere but is regular everywhere outside the stationary 
limit surface. We expect that, where the states discussed 
above are divergent, it is because the states fail to be 
Hadamard on that particular surface. However, our con- 
clusions are based on numerical computations only and 
we do not claim to have any rigorous results on the sin- 
gularity structure of the Green's functions defining the 
various states. 



2. Rate of rotation of the thermal distributions 

One of our key motivations for studying quantum 
fermion fields on Kerr space-time was to construct the 
analogue of a 'Hartle-Hawking' state, namely a thermal 
state. We have two candidates for this analogue state: 
our new state \H) (see Scc. lIII (J|) . and the state \CCH~) 
(see Sec. IIII A 3|) . These two states have some attractive 
regularity properties, as discussed in the previous sub- 
section. Given that the Kerr black hole is rotating, we 
now investigate the rate of rotation of the thermal distri- 
butions represented by the states \H) and \CCH~). 

To do this, we follow the method of [29j. Consider an 
observer moving on a world line with constant r and 9 
but with angular velocity 



= -f- 



d<p 
dt' 



We can associate a tetrad (e 



(t).e(r),e (e) ,e (¥j)i 



(4.34) 



with this 



observer. The vectors er r \ and erg\ are parallel to d/dr 
and 8/ 89 respectively, and the other two tetrad vectors 
are 1291: 



e (t) 



1_ (d_ 8_ 

M \dt + dip 



1 



1 



Af 



9t<p 

(gtt + togtip) 



9tt9ipip 

8_ 

dip 



(gt v + ty vv ) ^7 



d_ 

dt 



where 



9tptp I 



(4.35) 



(4.36) 



N = \9tt + 2£lgt<p ■ 

As well as the specific cases of a static observer (17 = 0) 
and a Rigidly Rotating Observer (RRO) with = Qh 
(|2.8[) , we are also interested in two non-constant values 
of Q. Firstly, if f2 = rizAMOj where 

9t<p 

3 



O 



ZAMO 



(4.37) 
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FIG. 8. The expectation values {T^u) CCH ~ B for components of the stress-energy tensor (multiplied by various powers of A 
(|2.2[) ; note that we do not claim that the power of A used necessarily corresponds to the rate of divergence of the components 
near the horizon), using the expressions ()C15h[C22|) with L = +1 (for L = —1, all components have the same values). The 
parameters used and format of the plots are the same as in Fig. 
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then the angular momentum of the stationary observer 
along the rotation axis of the black hole is zero. In com- 
mon with previous terminology Q , we call such observers 
Zero Angular Momentum Observers (ZAMOs). For com- 
parison with previous studies of the rate of rotation of a 
thermal distribution of spin-1 particles on Kerr J29[, we 
also consider a stationary observer with angular velocity 



n 



Carter 



r z + CL Z 



(4.38) 



whose orthonormal tetrad (|4.35[) is the Carter tetrad |56| . 

Following [29| , we study the angular velocity of an ob- 
server such that (T(t)( ¥ ))) = 0, where we are considering 
the expectation value of the stress-energy tensor operator 
in the state of interest. Such an observer is, in the termi- 
nology of j29| , a Zero Energy Flux Observer, who sees no 
angular flux of energy in that state. The angular velocity 
of a Zero Energy Flux Observer is denoted fizEFO- The 
angular velocity fizEFO can be computed from the ex- 
pectation values of the components of the stress-energy 
tensor in that particular state, as follows. The condition 
(T(t)(<p)) — means that fizEFO satisfies the following 
quadratic equation: 



where 12 £ 



An 2 ZEFO + Bn ZEFO + C = 0, (4.39) 

B = g<p<p(T tt ) - 9tt{T vv ), 

C = g tv (f tt )-g tt (f tv ). (4.40) 

In order to minimize numerical errors near the event hori- 
zon, the solution of this quadratic equation is written as 

M 



2C* 



ZEFO 



B± VB 2 -4AC 



(4.41) 



where the sign before the square root is chosen so that 
^zefo is regular and positive. 

In Fig. [§] we plot fizEFO for the expectation values 
(f^) CCH ~- B ~ (left-hand-plot) and (f^) H - B (right- 
hand-plot), together with f2# (|2.8[) (the angular speed 
of an RRO), fi ZAM o COB and Charter (j4^8| . From 
the left-hand-plot of Fig. U we see that (f^) CCH ~- B ~ 
is rigidly rotating close to the event horizon, but that 
the rate of rotation decreases as we move away from the 
horizon. Away from the horizon, the rate of rotation is 
slightly larger than both fizAMO and ^Carter- These re- 
sults are in qualitative agreement with those found in (29| 
for the electromagnetic case. It is the reduction in rota- 
tion rate as we move away from the event horizon which 
enables the expectation value (T^) CCH ~ B to remain 
regular everywhere outside the event horizon. 

The results for (T^ V ) H ~ B arc strikingly different. From 
Sec. IIV C 11 this expectation value is regular outside the 
event horizon and inside the speed-of-light surface. From 



Fig. [9] we see that, close to the event horizon, this expec- 
tation value is also rigidly rotating with the same angular 
speed as the event horizon. As we move away from the 
event horizon, rather surprisingly the rate of rotation of 
this expectation value increases above that of the event 
horizon, although it does not deviate away from flu by a 
large amount. The rate of rotation remains greater than 
fin until we reach the specd-of-light surface, where, from 
Sec. IIVC 11 the expectation value diverges. A stress- 
energy tensor which is isotropic and rotating rigidly with 
the same angular velocity as the event horizon is known 
to be divergent on the speed-of-light surface [57|, so the 
divergence of (T llv ) H ^ B on the speed-of-light surface is 
not surprising given that it seems to rotate a little quicker 
than Qh. 



V. PHYSICAL PROPERTIES OF THE STATES 

In this section we bring together our results and discuss 
the physical properties of the various quantum states we 
have defined. 



\B- 



: This state is defined in Sec. IIII A II as an absence 
of particles in the "in" modes at past null infinity 
I" and an absence of particles in the "up" modes 
at the past event horizon TiT . At future null infin- 
ity I + there is an outwards flux of particles in the 
super-radiant regime ujui < 0, corresponding to the 
'Unruh-Starobinskh' radiation [l5|,[l6[. The state 
is regular everywhere except on both the future and 
past horizons, where it diverges. It is not invariant 
under simultaneous t — ip reversal symmetry. 



\U~): To define this state (see Sec. IIII A 2[) . there are 
no particles in the "in" modes at X~ but the 
"up" modes are thermalized with respect to the 
frequency Co (corresponding to taking positive fre- 
quency modes with respect to an affine parame- 
ter along T-L~). This state is regular everywhere 
outside the event horizon. We expect that it will 
be regular on the future event horizon T-L + but di- 
vergent on the past event horizon H~ . Physically, 
this state corresponds to a star collapsing to form 
a black hole (for which space-time the past hori- 
zon H~ is unphysical so the divergence of the state 
there is not important). At future null infinity I + 
this state contains an outgoing flux of Hawking ra- 
diation. Like the 'past-Boulware' state \B~), the 
'past-Unruh' state \U~) is not invariant under si- 
multaneous t — ip reversal symmetry. 



\CCH~): This state is defined in Sec. IIII A 31 by adding 
a thermal flux of "in" particles, thermalized with 
respect to the frequency ui, to the 'past-Unruh' 
\U~ 



state 
\B~) and 
ncous t — 



271 ] . Like the other two 'past' states, 



\U ), it is not invariant under simulta- 
tp reversal symmetry. Due to this lack 
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FIG. 9. The rate of rotation fizEFO of the Zero Energy Flux Observer (ZEFO) (|4~4Tj) , for the expectation values (T^) 00 " ~ B 
(left) and (T M „) H_B (right). In both plots, we also show the angular velocity of the horizon Qh (|2.8[1 (denoted "rigid rotation"), 
and on the left-hand plot we also show Ozamo (|4.37p and ficartcr (|4.38[) . All quantities are plotted as functions of the co- 
ordinates (r,8). The expectation value in the right-hand-plot diverges on the speed-of-light surface, which can be seen in the 

numerical noise in the red surface. As in previous figures, we use the value a — cio = M <J 2 [a/2 — l] for the rotation parameter 
of the Kerr black hole. 



of time-reversal symmetry, the state \CCH~) can- 
not represent a black hole in a thermal equilibrium 
state, however it has a number of attractive prop- 
crtics, first noted in the bosonic case [HI, [29[. In 
particular, like \B~) and \U~), we have numerical 
evidence that \CCH~) is also regular everywhere 
outside the event horizons. We expect that it will 
be regular on at least the future event horizon Ji + 
as well. Close to the event horizon, the expecta- 
tion value (T^) CCH ~ B rotates with the same 
angular speed as the event horizon, but its angular 
speed then decreases as the distance from the event 
horizon increases. 

\B): This state is defined in Sec. IIIIBI bv an absence of 
"in" particles at past null infinity T~ and an ab- 
sence of "out" particles at future null infinity 1 + , 
which translates into an absence of both "in" and 
"up" particles far from the black hole. This state 
is therefore as empty as possible at infinity, and 
does not contain the outgoing 'Unruh-Starobinskh' 
radiation which is present in the 'past-Boulware' 
state \B~). However, the state \B) diverges inside 
the ergosphere. It is regular everywhere outside the 
stationary limit surface. Unlike the 'past-Boulwarc' 
state \B~), the state \B) is invariant under simulta- 
neous t — ip reversal symmetry. This is the natural 
vacuum state as seen by a static observer very far 
from the black hole. 

\H): This state is defined in Sec. MI CTl bv taking modes 
to have positive frequency with respect to affinc pa- 



rameters on the past and future horizons % . This 
corresponds to thermalizing both the "in" and "up" 
modes with respect to the frequency a). It is regular 
outside the event horizon up to the speed-of-light 
surface, where it diverges. We would anticipate 
that this state is also regular on both the future 
and past horizons Ti^. This state has some similar 
features to a rigidly rota ting thermal distribution 
of fermions in flat space |46[ which is also regular 
up to the speed-of-light surface. The state \H) is 
also invariant under simultaneous t— ip reversal. We 
conclude that our state \H) may represent a Kerr 
black hole in equilibrium with a thermal heat bath 
rigidly rotating with the same angular velocity as 
the event horizon. 

\B): This state, defined in Sec. MI 0"2l corresponds to an 
absence of "up" and "down" particles at the future 
and past event horizons Ti^. Like \H), it diverges 
on and outside the speed-of-light surface but is reg- 
ular inside the speed-of-light surface and outside 
the event horizon. We expect that it also diverges 
on the future and past event horizons H ± . Like 
both \B) and \H), it is invariant under simultane- 
ous t—ip reversal. Physically, this state represents a 
'rotating vacuum', that is, it is the vacuum state as 
seen by an observer rigidly rotating with the same 
angular velocity as the event horizon. This inter- 
pretation of the states \H) and \B) is borne out by 
the calculation of Ozefo in Sec. IIV C 21 where it 
is seen that the expectation value (T^) 11-3 corre- 
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sponds to a state which is rotating with almost the 
same angular speed as the event horizon. 



VI. DISCUSSION 

In this section we summarize the key results of this 
paper and discuss the wider implications of our work. 



A. Summary of our results 

In this paper we have studied in detail the quantum 
field theory of massless spin-1/2 particles propagating on 
a Kerr black hole. We began by reviewing the formalism 
for massless fermions on the Kerr geometry, and describ- 
ing the classical "in" and "up" field modes. The lack of 
super-radiance for fermionic fields, shown in Sec. Ill C[ is 
our first indication of a difference between the behaviour 
of bosonic and fermionic fields on rotating black hole 
space-times. 

In Sec. HIT] we tackle the subtle issue of quantizing the 
fcrmion field and constructing suitable quantum states, 
before numerically computing expectation values of the 
fcrmion current and stress-energy tensor for these states 
in Sec. IIVI In the absence of a methodology for calculat- 
ing renormalized expectation values on Kerr black holes, 
we have had to restrict our attention to differences in 
expectation values in two quantum states. 

We began with the uncontroversial 'past-Boulware' 
\B~) and 'past-Unruh' \U~) states which have been suc- 
cessfully constructed for bosonic fields. We also con- 
sidered the state \CCH~) [27j which is constructed by 
adding a thermal distribution of "in" particles to the 
\U~) state. All three 'past' states above can be defined 
for bosonic and fermionic fields, and all three are regular 
outside the event horizon. 

For bosonic fields, a 'Boulware' state empty at both 
future and past null infinity cannot be defined p8l |29| 
(see also App. IA 2[) . However, for fermionic fields we 
have been able to define such a state, \B). Unlike the 
'past-Boulware' state \B~), the state \B) is not regular 
everywhere outside the event horizon, but diverges inside 
the ergosphere. 

One of our original motivations for this study was 
the non-existence of a true 'Hartle-Hawking'-like state 
for bosonic fields on Kerr space-time [25j|. As well as 
the \CCH~) state discussed above, in the literature the 
state \FT) has been postulated to be an analogue of 
the 'Hartlc-Hawking' state for rotating black holes. For 
scalar fields, the state \FT) is regular only on the axis 
of rotation of the black hole [28[ . In this paper we have 
defined a state \H) which is the fermionic analogue of 
the state \FT). The state \H) is rather better behaved 
than the bosonic \FT) state, being regular between the 
horizon and the speed-of-light surface and divergent on 
and outside the speed-of-light surface. This state is the 
closest we have to a 'Hartle-Hawking' state for fermions 



on Kerr. Whereas Frolov and Thorne Q had to use an 
^-formalism to define their state, for fermions we are able 
to define the state directly by an appropriate definition 
of positive frequency. 

Finally, we have also defined a modified 'Boulwarc'-likc 
state, \B), which is empty as seen by a rigidly- rotating 
observer close to the event horizon. 

The regularity properties of all the states considered 
in this paper are summarized in Tab. HI and the physical 
properties of the various states are discussed in Sec. |V] 



B. The behaviour of bosonic and fermionic fields 
on Kerr 

A central theme in our work has been the differences 
between the quantum field theory of bosonic fields and 
the quantum field theory of fermionic fields on Kerr 
space-time, particularly in relation to the construction 
of 'Boulware' and 'Hartle-Hawking' states. At a clas- 
sical level, the fundamental difference between bosonic 
and fermionic fields is that bosonic fields exhibit the 
phenomenon of super-radiance for modes with frequency 
Oju) < (super-radiance means that an incident wave 
in this frequency range is reflected back to infinity 
with greater amplitude than it had initially). Super- 
radiance is a consequence of the weak-energy condition 
for bosonic fields. Classical fermionic fields do not obey 
the weak-energy condition and do not exhibit super- 
radiance, meaning that fermionic waves incident on a 
rotating black hole are always reflected back to infin- 
ity with an amplitude no greater than the incident am- 
plitude. At first sight it is not clear what the conse- 
quences of this classical phenomenon are for quantum 
field theory, particularly when its quantum analogue (the 
'Unruh-Starobinskii' effect [l5|, corresponding to sponta- 
neous emission in those modes with Cjlu < 0) occurs for 
both bosonic and fermionic fields. 

The proof of the Kay-Wald theorem [25[ on the non- 
existence of a 'Hartle-Hawking' state for scalar fields on 
Kerr uses an energy condition which arises from super- 
radiance. This may indicate that classical super-radiance 
plays a deeper role in the quantum field theory of scalars 
on Kerr, but there is no indication that it is a necessary 
condition for the non-existence of the 'Hartle-Hawking' 
state. Indeed, in this paper, for fermions we have not 
been able to construct a state satisfying all the conditions 
of the Kay-Wald theorem (regularity everywhere on and 
outside the event horizon and respecting all the symme- 
tries of the space-time) although fermions do not have 
classical super-radiance. Of course, this docs not mean 
that such a state does not exist, but the natural defini- 
tions do not yield such a state and we suspect that an 
analogue of the Kay-Wald theorem does hold for fermion 
fields (although proving such a statement would not be 
straightforward) . 

At a technical level, the existence of super-radiant 
modes appears to make the quantization of bosonic fields 
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at horizon 


inside ergoregion 


outside SoL 


IB'} 


X 


/ 


/ 


\u-) 


/ (on H + only) 


/ 


/ 


\CCH~) 


/ (on H + only?) 


/ 


/ 


\B) 


X 


X 


/ 


\H) 


/ 


/ 


X 


\B) 


X 


/ 


X 



TABLE I. Regularity properties of quantum states for fermions on a non-extremal Kerr black hole. A / indicates that the 
state is well-defined in this region, whereas a X indicates that it is divergent. The notation SoL means 'speed-of-light surface'. 



more complicated (see, for example, @, H, H HI ) ■ The 
key reason for these technical difficulties is the need, for 
bosonic fields, for modes designated to represent 'parti- 
cles' to have positive norm, while those for 'anti-particles' 
must have negative norm, so that the canonical commu- 
tation relations hold. For bosonic fields, this greatly re- 
stricts the possible choices of positive frequency used to 
define quantum states. In particular, in the terminology 
of Sec. Ill CI the "in" modes for bosonic fields have posi- 
tive norm if lo > 0, while the "up" modes have positive 
norm if u> > 0. This causes problems in defining quan- 
tum states (see related discussion in App. I A 2[) . For ex- 
ample, in Sec. IIII Bl we construct a candidate 'Boulware' 
vacuum \B) for fermions, which is empty of both "in" 
and "up" mode particles with frequency cj > 0. Such a 
construction is not immediately possible for bosons due 
to the need to have Q > for the "up" modes. This 
problem can be circumvented to some extent for bosonic 
fields by the use of, for example, the r/ formalism of Q, 
but, as discussed earlier, the resulting states have some 
unattractive features. 

The key difference between bosonic and fermionic 
fields, as far as the definition of quantum states is 
concerned, is that all fermionic modes have positive 
norm (due to the fact that fermionic fields satisfy anti- 
commutation rather than commutation relations). Wc 
arc therefore free to split the quantum field into pos- 
itive and negative frequency modes without worrying 
about the norm of those modes. This provides much 
greater freedom in the choice of quantum states, as seen 
in Sec. IIIIl Of course, it does not guarantee that any of 
those states are physically reasonable nor the regularity 
of those states, but the fact that there is so much more 
freedom in defining states for fermions compared with 
bosons means that one is more optimistic about being 
able to find states which have attractive physical proper- 
ties. 



C. Broader issues 

This paper has been concerned with the quantization 
of massless fermion fields on a non-extremal Kerr black 
hole. In this section we conclude our discussions with 
some initial thoughts on the application of our results 
to the alternative situations of a massive fermion field 



and/or an extremal Kerr black hole. 

1. Massive fields 

In one sense the inclusion of fermion field mass would 
represent a technical complication in our analysis (in 
particular, the upper and lower two-spinors in our four- 
spinor (|2 . 22[1 would no longer be proportional), but it 
could also change the underlying physics. Due to super- 
radiant scattering, massive scalar fields have unstable 
bound states with energy < u < mfijj, which pro- 
duces the 'black hole bomb' effect (see, for example, 
[59l - [6]| ). Since classical fermion fields do not exhibit 
super-radiance, a 'black hole bomb' effect is not antic- 
ipated for fermions. Instead, it has been suggested (62j] 
that massive fermion modes in the super-radiant regime 
(0 < lu < mf2jy) condense and form a so-called 'Fermi 
sea' surrounding the black hole. Although fermions 
are subject to the quantum analogue of classical super- 
radiance, namely the Unruh-Starobinskii fl5llT(l| sponta- 
neous emission of particles in the super-radiant regime, 
the Pauli exclusion principle means that there can be at 
most one fermion in each state, preventing the exponen- 
tial build-up in the 'black hole bomb' scenario. It should 
be emphasized that the 'black hole bomb' is a classi- 
cal effect, while the proposed Kcrr-Fermi sea would be 
primarily quantum in origin. It would be interesting to 
investigate in detail the quantum field theory of a mas- 
sive fermion field on a Kerr black hole, and elucidate the 
effect of the Kerr-Fermi sea on the quantum states we 
have defined in this paper. 

2. Extremal Kerr black holes 

Quantum fields on an extremal Kerr black hole are 
of central importance for the Kerr-CFT correspondence 
[H (see also [HS1 for reviews). The Kerr-CFT corre- 
spondence is concerned with the near-horizon geometry 
of an extremal Kerr black hole. The extreme Kerr ge- 
ometry has metric l|2.ip . with a = M so that there is a 
single (degenerate) horizon at r — M with zero Hawking 
temperature and horizon angular velocity 
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The near-horizon geometry is obtained as a scaling limit 
of the metric (|2.ip by defining new co-ordinates as follows 



in the extremal limit the Boltzmann factor (16.41) becomes 



t — >• t = y, r -> r = A/ + Xr, 
X 



Lf ^ If = If + 



t 



2MX 
,3.2) 

and then taking the (well-defined) limit A — >• 0. The 
resulting metric is no-longer asymptotically flat, but re- 
sembles AdS2 x S 2 . A particularly important feature of 
the near-horizon geometry for the Kerr-CFT correspon- 
dence is that it has an enhanced symmetry compared 
with the Kerr metric, namely a third Killing vector of 
the form 



Co = rdr - idj, 



(6.3) 



giving an SL(2,R) x U(l) isomctry group, which is ex- 
ploited in the CFT part of the correspondence. 

For the CFT correspondence to make sense, and, in 
particular, for the counting of the microscopic CFT states 
to yield the classical entropy of the extremal Kerr black 
hole, it is necessary for the CFT dual to have a non- 
zero temperature. This implies that there is a non-zero 
temperature for the state of a quantum field on the near- 
horizon geometry, which in turn requires a suitable defi- 
nition of a thermal state on the full extremal Kerr geom- 
etry. Even though the Hawking temperature of the ex- 
tremal Kerr black hole is zero, such a temperature is de- 
fined [HI through first considering a near-extremal Kerr 
black hole. 

Assuming for the moment that such a state can be 
constructed, suppose that a thermal state is defined for 
a near-extremal Kerr black hole with temperature Tjj . In 
this state quantum field modes are thermally populated 
with Boltzmann factor 1321 



cxp 



H 



(6.4) 



The question is then how to take the near-horizon ex- 
tremal limit. In order to see how this could be done, we 
need to consider the field modes in more detail. 

For a massless "up" field mode (see (|3.8p or (|Al|l). 
taking the near-horizon limit (|6.2[) corresponds to con- 
sidering only modes with ui = mflH = m/2M on the 
full extremal Kerr geometry (63|. These modes are 
rather special, lying on the boundary between the super- 
radiant and non-super-radiant regimes for bosonic fields 
(recall that fermionic fields do not exhibit classical super- 
radiance). Such modes are contained within the region 
close to the event horizon and are decoupled from the 
asymptotic region of the full extremal Kerr geometry [63j . 
This is consistent with reflecting boundary conditions at 
the AdS-likc boundary of the near-horizon geometry. 

The extremal limit of the Boltzmann factor (|6.4p can 
now be taken as follows. Setting oj = m/2M and defining 
the 'temperature' T v by 34 



T, 



v Th 

t™o (l/2M)-f2 H : 



(6.5) 
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T 



(6.6) 



The CFT interpretation of this temperature (and, in- 
deed, the taking of the extremal limit) do not concern us 
here. Rather, we comment on the sense in which a ther- 
mal state can be defined for a near-extremal Kerr black 
hole (which is central to the definition of the extremal 
temperature above). 

As recognized in the Kerr-CFT literature [H, [HI, the 
fact that the Kerr metric does not possess a globally 
time-like Killing vector makes defining thermal states 
rather difficult. Nonetheless, the non-extremal Boltz- 
mann factor (|6.4I) is justified in the Kerr-CFT literature 
as a 'Frolov-Thorne temperature' [H, HH coming from 
the Frolov-Thorne state \FT). As discussed elsewhere in 
detail (see [28[ and Add. |A"1 below), for bosonic fields the 
Frolov-Thorne state is regular only on the axis of sym- 
metry and therefore is ill-defined even on regions very 
close to the horizon. We have seen in this paper that the 
analogue of the Frolov-Thorne state for fermionic fields 
is regular outside the event horizon and inside the speed- 
of-light surface, so it may be possible to justify the CFT- 
tcmperature using a quantum state of thermal fcrmions 
close to the event horizon of a non-extremal black hole. 
However, we will not explore this possibility further in 
this paper. 

A related question is whether it is possible to define 
sensible quantum states directly on the extremal Kerr 
black hole geometry (either the full geometry or the near- 
horizon geometry). Given that the extremal Kerr black 
hole does not possess a globally time-like Killing vector 
(neither the near-horizon limit [34| nor the full space-time 
[3a|). one might anticipate that many of the challenges 
of defining quantum states on non-extremal Kerr black 
holes would remain. At first sight, it would seem that 
the fact that the Hawking temperature of an extremal 
Kerr black hole is zero might simplify matters. However, 
for extremal Kerr, the speed-of-light surface crosses the 
horizon at a latitude 9 = arcsin (\/3 — I) [13, [HI- It is 
therefore difficult to envisage how even a rotating vacuum 
state (the analogue of our \B) state) might be defined 
in the extremal case. We leave this interesting question 
open for future investigation. 
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Appendix A: Quantum field theory of scalar fields on Kerr 

In this Appendix, for ease of reference, we briefly outline some of the key features of scalar quantum field theory 
on Kerr. The notation follows (28j . where further details may be found. 



Scalar modes 



A complete orthonormal basis of mode solutions of the Klein-Gordon equation is defined, for to > 0, as follows 
Sill: 



1 



u A = 



up _ 



V8tt 2 o; (r 2 + a 2 ) 
1 

y/8ir 2 u (r 2 

1 

^8tt 2 (-Co) (r^+a 2 )' 



r iut e im vS A (9)R%(r), Co > -m£l H , 
i- iut e im vS A {9)Rf(r), u> > 0, 



a iuit —imip i 



S A (0)i£ p A (r), > Co > -mn H , 



(Al) 



r 



where A = {uo,i,m}, —A = {— — m}, the functions 
5a are the usual scalar spheroidal harmonics, and the ra- 
dial functions R l ^ np {r) have the asymptotic behaviours: 



then promote the expansion coefficients a a to operators: 



(A3) 



In order that the operators a A satisfy the canonical com- 
mutation relations 



pup iur* 



— oo 
oo 



r* — > — oo 



(A2) 



We remark that care has to be taken in the definition of 
the "up" modes in (|A1|) because of the need to consider 
only modes with positive norm. The norm of the "up" 
modes is proportional to Co, meaning that we have to 
consider separately those modes with Co > and Co < 0. 
This is one of the subtleties which plagues scalar quantum 
field theory on Kerr space-time. 

The reason for considering only positive norm modes 
is the following (45[. We wish to expand the quantum 
scalar field as a sum over modes (compare (|A9[) ) and 

I 



^AA' 



(A4) 



it must be the case that the modes u A have positive norm 
and the modes u A have negative norm. This restricts the 
way in which candidate vacuum states can be defined 
(see, for example, [45[ for the simpler case of rotating 
Minkowski space). 

The following relations between the coefficients in (|A2|> 
hold: 



I A in I 2 w I r>ir |2 

0A\ — — 0-Dy\ 
LO 



'A 



\oATf = U 



^ios A p i 2 : 
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, , R in * A U P — ,~< R U P A in * 
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Lo B i £ = oo Bl P . 



(A5) 



The first two of these relations show that for loCo < 0, 

both o-^a ancl lo^- A P | are greater than unity, indicat- 
ing super-radiance. 



An alternative, complete, orthonormal, set of basis modes can be defined for to > [2£ 

rrnf 1 



^87^0; (r 2 + a 2 ) 
1 

^87T 2 w (r 2 + a 2 ) 
1 

y/8TT 2 (-Cb) (r 2 +a 2 ) 



e -™t e ™<f>S K {9)R l Z*{r), Co > -mn H , 
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(A6) 
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Both the u° A nt and u A ° wn modes can be written in terms 



of the 



and u A p modes. 



For non-superradiant modes 



(w > 0,ui > 0), the results are: 



out 
u A 

down 
l A 



A in* in 

a A A U A 



W rjin* in 
— 0« A U A 
LU 



in* up 



up* up 
° l A > 



(A7) 



and it should be noticed that the right-hand-sides of these 
equations involve u A p and not their complex conju- 
gates. However, for super-radiant modes ujuj < 0, the 
situation is different: 



/(in* in 
O^A U A - 



i — j->up up* 
LU 



,down 




pin in* 



jup* up 



(A8) 



The important point about the relations (|A8[) is that 
they involve the complex conjugates of the "in" and "up" 
modes. This means that one obtains non-trivial Bogoli- 
ubov coefficients for super-radiant modes when changing 
from a basis of "in" and "up" modes to a basis of "out" 



and "down" modes. The result of this is that the vac- 
uum defined using the "in" and "up" modes as a basis 
(the 'past-Boulware' state \B~) defined below) is not the 
same as the vacuum defined using the "out" and "down" 
modes as a basis as far as the super-radiant modes are 
concerned. This is precisely the phenomenon of Unruh- 
Starobinskii radiation. 



2. Defining quantum states 



The 'past-Boulware' state \B ) is defined by first ex- 
panding the scalar field in terms of the u A and u A p basis 

(|A1|) and promoting the expansion coefficients a A ^ up to 
operators satisfying the usual commutation relations: 



E 
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1=0 m=-t 
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\ up~up 
[ U A a A 
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Then the 'past-Boulware' state is defined as the state 



annihilated by the operators d' 



l/up 



The definition of the 'past-Unruh' state \U~) could, in principle, follow that in Sec. IIII A| but the super-radiant 
modes, coupled with the need to use only positive norm modes (so that the "up" modes (|A1|) are only defined for 
Co > 0) complicates matters. We do not present a full derivation here, as it can be found in Appendix B of ML The 
simplest way to illustrate the nature of the resulting state is to give the expression for the two-point function 0, HH : 

G u -{x,x') = (U-\&{x)${x , )\U-) 

= E E I r^coth(-^-)u A p (x)uT(x')+ fd^KxXtx')}, (A10) 



from which it is clear that the "up" modes are thermally 



Now suppose that we wish to attempt to define a 'Boul- 
ware' state empty at both X~ and T + . Such a state 
would need to be constructed from the "in" and "out" 
modes (see (|A1|) and (|A6|) respectively) and it would be 
the boson equivalent of the fermion state \B) defined 
via Eq p.20[) . Such a state was suggested some time 
ago [6J], although its properties have not been investi- 
gated. The "in" and "out" modes are not orthogonal, 
and so we would need to write the "out" modes in terms 
of the "in" and "up" modes, using the relations (|A7HA8p . 
The resulting coefficients of the creation "up" operators 
then turn out to have positive norm in the superradiant 
regime (see Eq. (6.3.3) in [65|), and so they should in fact 
be annihilation operators. We could therefore make use 



populated. 

I 

of the 77- formalism introduced by Frolov and Thorne 0] ■ 
However, the FT-state (see Eq. (|A12I) below) constructed 
in Q using the 77-formalism is actually ill-defined every- 
where (except on the axis of symmetry). It is therefore 
likely that the 'Boulware'-like state that we have just 
suggested for bosons, even if formally empty at I~ and 
T + , is similarly ill-defined in most of the space-time; we 
leave such a question for future investigation. 

For scalar fields, the theorems of Kay and Wald [25| 
prove that there does not exist a Hadamard state on 
Kerr space-time which respects the symmetries of the 
space-time and is regular everywhere. In the absence of 
a 'true' Hartle-Hawking state as a consequence of this 
result, there have been a number of attempts in the lit- 
erature to define a 'Hartle-Hawking'-like state. 



The first such attempt is due to Candelas, Chrzanowski and Howard [27J, where the "in" and "up" modes are each 
thcrmalizcd with respect to their natural energy, so that the two-point function for a scalar field in this state is given 



31 



by 

G C ch-(x,x') = {CCH-\&(x)*(a/)\CCH- 



oo I 

EE 

£=0 m=-t 



r ^ coth (JL-^J u%(x)uT(x') + £ dCo coth (^j u n*K p v)} ■ ( Al1 ) 



It is argued (at least for scalar fields) in J28[ that the CCH-state \CCH~) is workable but does not represent an 
equilibrium state. In particular, it is not invariant under the symmetry transformation (t, ip) — > {—t,—ip) of the 
underlying Kerr space-time. Detailed calculations of the differences in expectation values of the stress-energy tensor 
for electromagnetic fields in the CCH-state and 'past-Boulware' state are presented in (29[. It is found that, close to 
the horizon, such differences correspond to minus a thermal distribution rigidly rotating with the event horizon, but 
that this rigid rotation does not seem to hold further away from the event horizon. No divergences in the CCH-state 
were found. We conclude that while the CCH-state has some interesting properties and appears to be well-behaved, 
it does not represent a black hole in equilibrium with a thermal bath of radiation at the Hawking temperature. 

A second candidate 'Hartle-Hawking' state was proposed by Frolov and Thorne @, and differs from the CCH-state 
in the choice of thermal factor for the "in" modes: 

G FT {x,x') = (FT\$(x)$(x')\FT) 

= EE {f~duco&(jlr)uAW<*W (A12) 



The FT-state \ FT) has the advantage over the CCH-state 
of being, at least formally, invariant under simultaneous 
t — ip reversal. However, it is argued in [28j that the FT- 
state is fundamentally flawed, and is regular only on the 
axis of rotation. Note that for scalars, one cannot replace 
the integral over u in (|A12[) with an integral over Cu be- 
cause the "in" modes are defined for u > 0, not Q > 0. 
Therefore we cannot, for scalars, define a direct analogue 
of the state \H) defined in Sec. IIII Cl for fermions. 



Appendix B: Dirac and spinor connection matrices 

In this Appendix we list the Dirac and spinor connec- 
tion matrices for the Kerr geometry using our space-time 
conventions. 



Dirac matrices 



where the flat-space 7° matrices are given by 



7° 



(ih \ 









I -ih ) ' 








-iaj 
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(B2) 



with I2 the 2x2 identity matrix and <Ji the usual 2x2 
Pauli matrices 




As anticipated, the flat-space 7° matrices (|B2j) satisfy 

{f\7 fc } =2if b . (B4) 
We also define a chirality matrix 7 5 by 



(B5) 




A suitable basis of 7^ matrices for the Kerr metric 
(|2~Tj) can be found in [lj|: 
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■7° 
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7° + 
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:T 2 , 
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2. Spinor connection matrices 

The spinor affinc connection matrices T M are most eas- 
ily computed by using a vierbein such that 

I" 1 = (B6) 

where 7° are the flat-space Dirac matrices (|B2|) . In terms 
of vierbein c omp onents, the spinor connection matrices 
are given by (23, E3, El] : 



(B7) 
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Using this formula we find that the spin connection ma- 
trices can be compactly written as follows [13] : 



r, 
rv 

r s 
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r 2 - a 2 cos 2 f I 7"7 
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(B8) 



where 



y4 = AS + 2Mr (r 2 



B = a 2 r cos 2 i 



a 2 M cos 2 9 + r + Mr 2 



(B9) 



These T v matrices also satisfy the additional condition 
Trr„ = 0@. 



Appendix C: Stress-energy tensor components 

The classical stress-energy tensor (|2.19p for a fcrmion mode tf>\ (here we omit the superscripts "in/up" because the 
formulae apply equally well to all modes) is: 

y 

aT^ = - [^a7mV„Va + V , a7^V m Va - (V^V'a) iAa - (V„Va) I^a] ■ (CI) 

In analogy with the quantity (|4.8p defined for the number current, we define the following quantity At^v , which is 
required for the computation of expectation values: 

Atytv — —aT^v — aT^ v . (C2) 

The expressions for the components of sT^v and At^u are rather lengthy and given below, where, for conciseness, 
we omit the subscript a and also all "in/up" mode labels. The notation 5ft denotes the real part and 9 denotes the 
imaginary part of complex functions. We have explicitly verified that these stress-energy tensor components satisfy 
the conservation equations V^T^^ = 0. The conservation equations for a classical stress-energy tensor on a Kerr 
space-time can be found in 1281 . although we note that there is an error in one of their equations. The v = t, 9 and 
f conservation equations in [28| are correct, but the v = r equation should read: 

d r (S T/) + — — d e (EsinflT/) - rT e e - A" 1 (ra 2 sin 2 9 - T) T r r 

= — [-T T u + 2aT sin 2 9 T tv + sin 2 9 (-Ta 2 sin 2 9 + rS 2 ) T w ] , (C3) 

where 

T = M (r 2 - a 2 cos 2 9) . (C4) 

The expressions (|C5HC24|) given below depend explicitly on L. The differential equations (|2.25|) satisfied by the radial 
functions also depend on L. The boundary conditions on the radial functions for the "in" (|2.32|) and "up" (|2.33j) modes 
are stated for L = +1 only. For L = +1 therefore, the radial functions satisfying the appropriate boundary conditions 
can be substituted into the stress-energy tensor components (|C5HC24|) . For L = — 1, the simplest way to obtain the 
corresponding expression for the stress-energy tensor components is to substitute L = — 1 into (|C5I IC24| and make 
the swap iR\ O 2-Rai since the differential equations (|2.25[) satisfied by the functions iR\ and 2-Ra swap over under 
the map L —L. The radial functions for L = +1, satisfying the original, L = +1, boundary conditions (|2.32[|2"T3"3"1) . 
can then be used in the computation of the stress-energy tensor components. Indeed, it is straightforward to see that 
the quantities (jC15flC24|) below used in Sec. IIVI in the computation of expectation values of the stress-energy tensor 
are invariant under the map L — > — L. 
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Firstly, we give the expressions for aT M!/ : 
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Secondly, we give the expressions for a£ m „, derived from those for aT^v using the symmetries (|2.28l[2~29l) : 
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+4a£L sin(9 [3? { t R* 2 R') - 3? Qi?' 2 R*)] X S 2 S 
La cos 



2Y?L 



A 



|ii?| 2 - | 2 i?| 2 [^ 2 + 2 5 2 ] 



r 2 +a 2 (l + sin 2 0)] \ X R\ 2 - | 2 i?| 2 [^ 2 - 2 S 2 ] 



|iii| 2 -| 2 /i| 2 [iS^-aS 2 ], 



/A 
" 8tt 2 £ 2 

— 2 J- v3 . - {-2VAS 2 (a W sin 2 + m) l^ 2 + \ 2 R\ 2 } ^S 2 + 2 S 2 } 
lD7r^v Az j j sin L J 

-8L£ 2 [(r 2 +« 2 )w + dm] sinf? S (ji? 2 iT) i5 2 5 

-VAcos6» [S 2 + 4Aira 2 sin 2 6»] \\ X R\ 2 + \ 2 R\ 2 ] [ X S 2 - 2 S 2 ] 

-4sin6» [(r - M) E 2 - 2M (r 2 + a 2 ) (r 2 - a 2 cos 2 6)} 5R (iR 2 R*) iS 2 S} , 

r \2ar sin 6 $t( 1 R 2 R*) + y/ASL [3 - 3 ( 2 iT 2 i?')] [iS" 2 + 2 S' 2 ] 

47r 2 A2Ssin6' 

-iaTAcose [|ii?| 2 + | 2 i?| 2 ] [i5 2 - 2 S 2 ] 



1 



8vr 2 VAS] sin 
1 



-4acos6»3?(ii? 2 i?*) xS 2 S ■ 



ar sin ( 



|ii?| 2 + | 2 i?| 2 [ 2 S 2 -iS 2 ] 



16tt 2 VAE 2 sin 
2E 2 

mL 



{2av / AS sin 2 [3 ^R* ji?') + 3 ( 2 iT 2 fl')] [i^ 2 + 2^ 2 ] 
ii?| 2 - | 2 i?| 2 l [i^ 2 + 2 S 2 ] - 4L (r 2 + a 2 ) Esin0 [3? (ii?* 2 i?') - 3? (i# 2 i?*)] iS 2 S 



L 
A 



a 2 A sin 2 
1 



(r 2 + a 2 ) 2 + 2M ra 2 sin 2 ' 



cos 8 



\iR\ 2 -\ 2 R\ 2 [ 2 S 2 -iS 2 ] , 



{2EL [%( 1 R 2 R*)-%( 1 R* 2 R)} ^S' 2 S - i5 2 S"] - 4rasm8$t {iR 2 R*) i5 2 5 



8tt 2 VAE sin 
-av^cosfl [|ii?| 2 + | 2 i?| 2 ] [i5 2 - 2 5 2 ] } 

1 

"47r 2 VAE 3 



AS m sin ( 



|ii?| 2 + | 2 i?| 2 



US* +2^1 - 4^mL (r* + a 2 ) 3 (ii? 2 i?*) iS.S 



ra sin t/ cos 

2 /a IY„2 „2 



lli?l 2 + l 2 i?l 2 [ls 2 - 2 s 2 ] 



-2Mo sin 2 6 [(r 2 - a 2 ) £ + 2r 2 (r 2 + a 2 )] 3? {iR 2 R*) iS 2 S) . 
These quantities arc used in the numerical computations in Sec. IIV CI 
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